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ABSTRACT 

This work reviews the existing literature on slam force 
prediction methods and discusses some of the theory of ship 
bottom slamming. One of these prediction methods is chosen 
as the best and further improvements are made on it by using 
a new model of slamming. 

The work is designed to determine a slam forcing function 
suitable for studying ship hull girder dynamic response. 

Slam forces are determined by integrating the pressure time 
history over the portion of the hull subject to slamming. 

A computer program is included which determines the slam 
forces and impulses from the ship's lines, principal dimen- 
sions, and seakeeping data. The method is applied to two 
ships, the MARINER and the FOTINI-L, and the results 
presented. 
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I. 



INTRODUCTION 



The primary objective of this thesis is to review the 
existing literature on slamming and to develop a method of 
predicting slam forces to determine hull girder response. 

This project is a portion of an overall study presently under- 
way at M.I.T. sponsored by the American Bureau of Shipping 
(A.B.S.) on hull stiffness criteria. 

Presently, explicit criteria for ship stiffness do not 
exist. The ship classification societies have controlled stiff- 
ness or deflections by limiting allowable Length to Depth (L/D) 
ratios [1]*. With the introduction of new materials such as 
high strength steels and aluminum in ship construction, the 
question of stiffness and what should be allowed becomes more 
important. 

Considering the present L/D ratios as hull girder stiff- 
ness criteria, what other factors that affect stiffness could 
be used in design? Bending stiffness of a ship's hull can be 
varied within reasonable limits by the following means: 

1. By changing the basic hull depth, 

2. By designing to higher stress levels with higher 
strength steels, and 

3. By using materials with different elastic moduli. 

Since other factors can more properly be used to control stiff- 
ness other than L/D, it is not satisfactory to use this as a 

*Numbers in square brackets designate References at the end 
of the paper. 
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criterion for design. Probably a new criterion based 
specifically on a measure of stiffness should be used. 

Before any quantitative information could be developed 
for limiting hull stiffness, factors which would be changed 
by stiffness variations were assessed [2] . Assuming strength 
criteria are maintained, stiffness variations will result in 
changing factors which determine hull girder dynamic response. 
In the Reference [2] report it was concluded that the factors 
of ship stiffness which should be of concern were; 

1. Slamming response for its contribution of stress 
components, especially amidships, 

2. springing for its cumulative stress effects, and 

3. steady-state propeller-excited vibratory motions 
and their deleterious effects upon personnel and 
main machinery components. 

Although propeller-excited vibratory motions will change with 
hull stiffness variations, it is felt that these problems can 
be dealt with by more efficient and less costly methods than 
changing the hull girder stiffness. Therefore the primary 
factors of concern are slamming and springing response with 
variations in hull stiffness. 

In the overall study the Kline-Clough lumped-mass vibra- 
tion program, S.H.V.R.S., has been used to analyze the ship 
response to slamming and springing [3] . Several ships have 
been looked at in great detail. Some of the results have been 
reported in references [3] and [4]. To date, in this work on 
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slamming, a constant slam impulse of 100 Ton-seconds has been 
used in the response analysis modeled as a half sine pulse 
applied at a single location. Results, such as stress levels, 
other than trends in response due to changes in stiffness, 
could not be determined until a better model of the slam force 
and a satisfactory means of evaluating force magnitudes from 
ship to ship could be found. 

The purpose of this work is to develop a model of the 
bottom impact slam force compatible with the S.H.V.R.S. program 
which will require a minimum of manhours and computer time yet 
give adequate results. The bottom impact slam force developed 
here is a function of hull form, ship speed, and sea state. 

A computer program is developed which predicts the most 
probable slam force as well as an extreme force to be used for 
design. Other slam information is also available such as 
frequency of slamming and slam pressures useful in designing 
hull plating. 
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II. BOTTOM IMPACT SLAMMING 



For this study the bottom impact slam is defined as the 
"heavy blow" resulting from the forebody hitting the water 
surface at a high relative velocity. This is in contrast to 
the phenomenon of bow flare slamming which is a consequence 
of emersion of the bow flare and the resultant hydrodynamic 
and hydrostatic forces. Much of the work to date has been in 
the form of predicting the local pressures developed to study 
plating damage. That can be considered the micro-problem. In 
this work we are concerned with the macro-problem which in- 
volves hull vibrations and possibly large midship stresses. 

Oakley [5] has analyzed the micro-problem in great detail. 
There are many factors that have been observed to measurably 
affect the pressure history. These factors include the 
following: shape, velocity, and phasing of the surface waves; 

geometry of the hull and hull reaction characteristics; air 
and water compressibility; two and three-dimensional air and 
water flow; bubbles, spray, and even such things as marine 
growth. Due to the complexity of the problem and the random- 
ness of the sea, model testing and statistics play an important 
role in the solution to the slamming problem. 

Tick [6] developed the equation for the frequency of 
slamming which is dependent on the following: 

1. forefoot emergence statistics, 

2. relative velocity statistics, and 

3. the angle between the keel line and the wave 
surface at the instant of impact. 
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In most work to date the equation is simplified to remove the 
angle dependence and only the requirements of (1) forefoot 
emergence and (2) a relative velocity greater than a threshold 
velocity are considered necessary for slamming. Neglecting 
this angular dependence *esults in overpredicting the frequency 
of slam impact as well as the slam force. A method of dealing 
with the angular dependence is presented in Section IV of this 
work. The two requirements for slamming are met by the ship 
in Figure 1. Using Tick's equation without angular dependence 
allows prediction of the frequency of slamming and is the basis 
of slam calculations in seakeeping programs [7] . The one 
necessary input to seakeeping programs is the threshold velo- 
city. How this is determined and the value to be used is 
discussed in Section IV. 

From Tick's work it is possible to conclude whether or 
not a ship is subject to slamming under a particular set of 
operating conditions, in particular the ship speed, sea state 
and draft at the station under consideration. Given that 
slamming does occur, the next step is to determine the slam 
force. Two basic methods have been used in the past: (1) 

integration of the slam pressure time histories over the slam 
area; and (2) determine an analytical solution based on the 
rate of change of momentum. In Section III both methods are 
discussed and one selected for this study. 



15 



I 




I 

I 









I 





Before comparing the two methods, some discussion is in 
order with respect to the mechanism of slamming. Oakley [5] 
considered the time history of slamming to be made up of 
three phases; 

1. Phase (I) is the time in which the body is 
approaching the free surface up until it makes 
contact. In this phase the air flow and the 
surface wave deflection are of predominant 
importance . 

2. Phase (II) is tha period when the body goes 
through the free surface and involves the problem 
of "contact" and more or less complete wetting of 
the hull. 

3. Phase (III) is the fully wetted problem in which 
the rate of change of added mass and the local 
surface deflection are dominant. 

A very flat body or "U" form hull may accelerate the free 
surface significantly during phase (I) while a more "V" 
shaped section may develop greatest forces only during phase 
(III). Figure 2 shows the relationship between the three 
phases. Much of the work to date has dealt with only one or 
two of these phases at a time. This is especially true with 
the many two and three-dimensional drop tests performed. 
Chuang [8] and others have demonstrated the important role 
of the air in accelerating the free surface in phase (I) and 
as a cushioning mechanism in phase (II) . The wedge-entry 
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problem has been studied by many authors, most recently by 
Chuang [9] , and agreement for large deadrise angles (phase 
(III) only) is good. The theory breaks down for low deadrise 
angles where phase (I) and phase (II) become important. 

Therefore any solution to the slam forcing function 
should take into account the forces developed during all 
three of the phases unless one phase can be shown to be much 
greater than the others. At this time it appears that neglect- 
ing phase (I) and phase ill) is only valid for "V" form sections. 
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III. SLAM FORCE PREDICTION METHODS 



During the last twenty years a great deal has been written 
on slamming, attesting to the difficulty of a satisfactory 
solution to the problem. Much of the work has been done along 
with ship vibration analysis techniques. This study does not 
go into the question of structural response. It is felt that 
adequate methods exist for analyzing the structural response 
of a ship if the various forcing functions can be adequately 
represented. The S.H.V.R.S. program described in Reference [3] 
which has been used in this overall study has been shown to 
give satisfactory results. Therefore the model of slam force 
was chosen with the S.H.V.R.S. program in mind. 

The slam force prediction methods are of two types: 

1. The slam force is considered as resulting from 
the rate of change of momentum of the added 
mass of the ship with a buoyancy correction 
term. This method will be called the momentum 
method . 

2. The second method consists of determining the slam 
pressure and integrating spatially as well as 
temporally over the slam area. This method will 
be designated the slam pressure method. 

Both of these methods were considered and the second was 
considered more valid and more easily applicable to the 
overall stiffness study. 
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RATE OF CHANGE OF MOMENTE METHOD 



As early as 1956 Szebehely [10] proposed that the slam 
force was a result of the rate of change of momentum of the 
added mass in the region of the slam. Further studies by 
Leibowitz [11] used this model of slamming to analyze the 
slam forces on a Dutch destroyer. Good results were obtained 
between the theory and full scale results for bow flare slams. 
More recently Mansour anc d'Oliviera [12], [13] have utilized 

the same slam force function to predict slam bending moments 
in regular waves. 

The slam force (P) is considered as the sum of an inertial 
term ^ buoyancy term (P 2 ) is given by 



P(x,t) = P^(x,t) + p 2 (x,t) 



D / N 
Dt<"^^r^ 



+ pgA 



at 



(mw^) + + PgA 



( 1 ) 



where m = m(x,t) is the added mass per unit length, w^ is the 
relative vertical velocity between the wave and the ship, and 
A is the sectional area. The total derivative is used (^) 
since the rate of change of momentum occurs not only temporally 
but also spatially. 

It should be noted that this is the phase (III) problem 
described by Oakley. No account is made for the force which 
comes during phase (I) and phase (II) of the slam process. It 
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is expected that this method should show good results for "V" 
form hulls and bow flare slamming only. The major drawback 
to this approach, even if it were valid for the bottom impact 
slam case, is it requires a large amount of seakeeping infor- 
mation to determine m(x,t). In reality the added mass must 
be determined as a function of draft at each station and 
frequency of encounter f c r each wave considered. To obtain 
this information requires a large expenditure in manhours and 
computer time. The end i esult is also valid only for regular 
waves . 

SLAM PRESSURE PREDICTION METHOD 

A great deal of the initial work in this area has been 
done to solve the micro-j roblem — that is, to determine the 
slam pressure locally foi plate design information. Much of 
the work in this area hai been done by M.D. Ochi, K.M. Ochi, 
and S.L. Chuang at the Ns val Ship Research and Development 
Center. A summary of th( slamming studies to date appears in 
Reference [14] . In this work by M.K. Ochi and Motter they 
develop the basic theory behind the approach to slamming pre- 
diction in this method. 

In this approach the pressure is considered as a function 
of velocity of impact and section shape only. The slam 
pressure is given by 

p = kf^ (2) 
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where k is a parameter to be determined from section form 
and f is the velocity of Impact. 

From the many two and three-dimensional drop tests as 
well as full scale measurements and model tests it has been 
found that equation (2) h elds quite well for slam impact pro- 
blems. Therefore once the form coefficient (k) can be deter- 
mined for each section, t le slam impact pressure can be 
determined if the velocit/ of impact is known. Since relative 
velocities are easily determined from seakeeping programs, 
knowing the value of k for each section allows for the deter- 
mination of the slam pressure which can be integrated over the 
hull form to determine the slam force. 

This second approach has the advantage that the form 
coefficient k is obtained from model data which includes the 
contribution from phase (I), (II), and (III) forces. In 

addition, recent work by Loukakis and Chryssostomidis [15] 
make available the seakeeping information necessary for slam 
studies for ships that can be approximated by Series 60 hull 
forms . 

The slam pressure approach was chosen since it requires 
a minimum of seakeeping information, is applicable to all 
hull forms, and the results are applicable in irregular seas. 
Unlike most engineering problems it also turns out that this 
approach is easiest to apply. 
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IV. PRINCIPLES AND PROCI'DURE OF THE SLAM PRESSURE METHOD 



Ochi and Motter ' s work described in Reference [14] pre- 
sents the theory and assumptions of the slam pressure method. 
In trying to apply the method it appeared that some of their 
assumptions resulted in overpredicting the slam force and im- 
pulse. In this section the theory and assumptions made by 
Ochi and Motter are discussed and modified where the assump- 
tions were not felt full^ justified. The following topics 
will be discussed in developing the method: 

1. Frequency of slam impact, 

2. Slam pressure ard hull form coefficient, 

3. Most probable extreme pressure, 

4. Extreme pressure for design, 

5. Sea state data end approximation for storm duration, 

6. Slam impact area of the hull, and 

7. Slam impact force and impulse. 

FREQUENCY OF SLAM IMPACT 

Tick [6] derived the equation for the frequency of slam 
impact at the bow which is dependent on the relative velocity, 
relative motion, and angie between the keel and the wave sur- 
face. By integrating over all angles to neglect the angular 
dependence. 
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where 



N = Number of slams per second 
s 



T = Ship draft at station under consideration 
f* = Threshold velocity 

= Twice variance of relative motion 
= Twice variance of relative velocity 
The number of impacts in a storm of duration t hours, N(t), 
is then given by 



Different values of the threshold velocity below which 
slamming does not occur have been used, but following Ochi ' s 
work [14] , 



will be used where = 12 feet/sec for a 520 FT ship and the 
Froude scaling law is applied. Because some of the slams 
predicted by equation (4) will occur above the angle which will 
cause slamming, equation (4) will overpredict the number of 
slams. This nimber is used to determine the extent of slam 
area and slam pressure. Overestimating N(t) will yield a 
larger slam area as well as a slightly higher pressure. The 
assumption is made that slamming will still occur forward of 
the area that is overestimated and that the pressure can be 
corrected by a relative angle correction. This will be 
discussed later. 



N(t) 




(4) 




(5) 
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SLAMMING PRESSURE AND HULL FORM COEFFICIENT (k) 



To estimate the slam pressure it is necessary to find a 
functional relationship for pressure vs. velocity. As mentioned 
earlier, it has been show i that equation (2) 

u *2 

p = k r 

holds very well and can b 2 used if the form coefficient, k, 
can be determined. 

Ochi and Motter [16] used mapping and regression analysis 
techniques to find the k 'alues from experimental data. This 
method was used to develo ? a computer program reported in 
Reference [16] to determi le the coefficients. This program 
was used as a subroutine .n this work but it did not give 
satisfactory results for ?road beam sections. The regression 
analysis had been perform ?d over relatively narrow ship sec- 
tions and the confidence "ange for the coefficients was 0.073 
to 0.190 with 10% risk of error. In applying their program 
two problems developed: 

1. The k values for broad beam flat sections were 
well outside the 0.190 confidence figure; and 

2. The mapping blew’ up in many of the broad beam 
cases . 

Since many ships of comme rcial interest could not be covered 
by the program, another m ithod had to be found for these wide 
forms . 
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M.D. Ochi [17] in work which was an attempt to correlate 
the k values from two and three-dimensional drop tests as well 
as model tests proposed that the following holds: 

b'2 

k = 0.02 ^ (6) 

where B' = half-breadth at 1/10 design draft 

A’ = half-section area at the 1/10 design draft 
This work, which was done prior to the regression analysis 
technique in Reference [16], was compared to the results ob- 
tained by the computer program. Better agreement inside the 
confidence range was obtained if a slightly higher constant 
was used. Therefore it is assumed that 

r-2 

k = 0.027 (6a) 

In the program developed in this thesis, equation (6a) is 
used when section values above 0.190 are obtained from the 
regression subroutine. This is done since equation (6a) yields 
larger k values and high pressures which is more conservative. 

Now that the form coefficients can be calculated, it is 
possible to determine the maximum pressure at each station 
along the ship if the velocity of impact can be found. 
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MOST PROBABLE EXTREME PRESSURE 



The impact relative velocity is approximated by the 
amplitude of relative velocity which is a random variable 
following the Rayleigh probability distribution [14] . Using 
equation (2) and the Rayleigh distribution function, the 
slamming probability distribution function is given by 

f(p) = X exp [-X (p - p^) ] (7) 



< 




where p = slamming pressure 



p = threshold pressure = kf 



.2 



X = l/(k Rp 

By applying order statistics to the probability distri- 
bution (7), the distribution of extreme value of pressure in 
n finite mamber of impacts, p^, is given by 



f(p^) = [n f(p) (F(p)} 



n 



= nXexp [-X(p^ - Pq)] {1-exp [-X(p^ - p^) ] ^ (8) 



p - p < «> 

'^o n 



where 

F(p) 



P 

/ f(p) dp 

Po 
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From equation (8) it is possible to derive the most 
probable extreme pressure in n-observations , p^. Figure 3 
shows the probability distribution function given by equation 
(8). The most probable extreme pressure, p^, can be obtained 
by setting the first derivative with respect to p^ equal to 
zero . 



d f(p ) 

— = 0 

d p 
' n 

Using equations (8) and (9), p^ is obtained 

p = k (f * + R! In (n) ) 

^n * r 



(9) 



( 10 ) 



This most probable pressure can be used to calculate slam 
forces and is one of two pressure options in the computer 
program . 

EXTREME PRESSURE FOR DESIGN 

Although it may be reasonable to compare p^ with the 
largest value observed in experiments, the probability that 
a pressure is higher than p^ is 0.632. Since the probability 
distribution function, equation (8) as plotted in Figure 3, 
is never zero at high pressures, the extreme pressure for 
design can be determined only if a small probability of being 
exceeded (a) is allowable. The extreme pressure for design, 
(a) , can be found from 
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FJGiJRE 3 - EXTREME PRESSURE PROBABILITY DENSITY FUNCTION 



I 







00 



/ 




( 11 ) 



p 



n 



Using equations (8) and (11) 



p^(a) = k [f^ - : n{l- (1-a) 



( 12 ) 



The extreme pressure for design is controlled by the pre- 
assigned small probabili -y, a , which can be specified by the 
designer. In the calculations in the computer program Ochi's 
[14] value of a = 0.01 is used, but can be changed if desired. 

SEA STATE INFORMATION 

In the calculations necessary to determine slam pressures 
the sea state must be specified to determine 

1. the relative motions and velocities and 

2. the duration of the storm 

in order to calculate the number of slams from equation (4). 

In order to make the application of the program as simple 
as possible, the only pa' ameter which needs to be specified is 
the significant wave hei< ht, H(l/3). This is made possible 
by the work done by Louk kis and Chryssostomidis [15] from 
which the r.m.s. relativ motions and r.m.s. relative velocities 
may be determined for tha forward stations using only the 
parameters of: 

1. C„ = Block coef icient 

2 . L/B = Length to Beam ratio 
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3. B/T = Beam to Draft ratio 

4. H(l/3)/L = S = Non-dimensional sea state 

5. F = Froude number = SHIP SPEED/ /gL 

These parameters are sufficient to give good initial design 
values for relative motion and velocity. There is also a 
dependence on longitudinal center of flotation and radius of 
gyration but it is not particularly sensitive to these 
parameters [15]. These tables were developed from Series-60 
hull forms and, therefore, can only be used for ships which 
can be approximated by thsse forms. This allows calculation 
of slam forces for these ships without running seakeeping 
programs. The option is ilso available to specify the relative 
motions and velocities if desired. 

The storm duration also needs to be known. Following 
Ochi and Motter [14], wave records have been analyzed in the 
North Atlantic. A plot cf the envelope of storm duration vs. 
significant wave height is shown in Figure 4. This envelope 
is used to determine the maximum duration to be expected for 
a ship operating at the significant wave height. This envelope 
is the default option in the program and is considered the ex- 
treme storm for design. If desired, a storm duration may be 
specified in the input d< ta to the program 
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FIGURE 4 - STORM DURATION VS SIGNIFICANT WAVE HEIGHT 
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SLAM IMPACT REGION 



Using the time of st orm fiom the previous section as well 
as the tabulated r.m.s. potions and velocities, equation (4) 
can be evaluated at each station considered to determine the 
number of slams during ar extreme storm. The number of slams 
is greatest at the bow ar i decreases rapidly as the midships 
crossection is approachec . The sections where N(t) - 1 is 
assumed to be the longiti dinal extent of slamming. 

From model tests anc drop tests [17] it has been found 
that the slam pressure it essentially dissipated at a depth 
of 1/10 of the design dr^ft, T' . Therefore the pressure needs 
to be integrated only ove r the forward portion of the hull 
where N(t) - 1 and below 1/10 design draft. 

SLAM FORCE AND IMPULSE 

The information is r ow available to determine the most 
probable and extreme pref sure for design p^ and the 

keel for each station as well as the slam impact area for a 
particular ship, speed, ^ nd sea state. To integrate temporally 
and spatially, more infoi mation is needed about the spatial 
variation of pressure in the girth direction as well as the 
time history of the slam. 

The local pressure ■( ime history as used by Ochi [14] is 
shown in Figure 5. The ( uration of the slam pressure, t , is 
given by 

= 0.1 /L/520 = 4.4 X 10“^ /L (13) 
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FIGURE 5 - SLAM PRESSURE TIME HISTORY 
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which is an average for a 520 FT vessel whether it is in the 
flat bottom region or the bilge. Froude ' s scaling law is used 
for ships of different leigths. The slam pulse is assumed 
triangular with equal ris i and fall times. The duration of 
slam is one area where it is felt additional work needs to be 
done. In Appendix A an a .ternate equation for t^ is considered 
and presented here as equ ition (14). 

tj^ = 3.6 X 10 ^ B' (seconds) (14) 

where B' = Half-beam a 1/10 design draft in feet 
Until further work can be done to experimentally determine 
which time duration is be ;t, either can be used. The computer 
program uses equation (13 but it is felt by this author that 
from the analysis of drop test results beam should be used as 
the scaling parameter ins .ead of length. 

The work done by Och assumes that there is no dependence 
on angle between the keel and wave surface, and that the worst 
slam case is when the sla » pressure travels the full extent 
forward or aft at the slowest traveling velocity. It is felt 
that the wave slope has a large effect in reducing the slam 
forces predicted by Ochi ' s slam model. Therefore a new model 
of slam is proposed which accounts for the dependence on wave 
slope . 

Numerous experiments have shown that the slamming becomes 
maximum in regular waves 'hen the ship encounters waves of 
length L on the order of 1.0 to 1.1 times the ship length, L, 
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[12], [18]. It is assume! that a similar condition has been 

approximated in irregular seas leading to the worst slam 
condition. For the relat .ve velocities to be maximum, the 
ship will be pitching dow i in the vicinity of zero pitch 
angle. For this to occur the slam will take place in the 
trough of the wave. Figure 6 depicts this model of slamming. 
The slam pressure will travel from the forward perpendicular 
aft and from the aft slam station forward meeting at what is 
called here the slam center. Assuming the wave encountered 
has height equal to the significant wave height allows calcula- 
tion of the slope between the keel and the wave surface. This 
angular dependence will be used to modify the slam pressures 
from what would be expected for flat impact. 

Chuang [9] performed a series of drop tests with varying 
deadrise angles. It is proposed here that the reduction in 
pressure with deadrise arjle can be used for similar angles 
in the longitudinal direction. Figure 7 is taken from 
Reference [9] and shows a significant reduction in pressure 
with deadrise angle. The region between 0 and 3 degrees is 
where the compressibility of the air comes into effect. It 
is assumed this is not as pronounced in the longitudinal 
direction due to the escape of air in the transverse 
direction. The higher ar gles are almost linear when plotted 
on paper and the slope remains almost constant for all 

impact velocities. From the linear region of Figure 7 a 
correction equation for pressure was found as 
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FIGURE 7 - MAXIMUM SLAM PRESSURE VS. DEADRISE ANGLE 
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(15) 



p^(a,0) = p^(a) 

where p^(a,6) = maximum keel pressure with angular correction 

Pn(“) = pressure predicted by equation (12) 

6 = angle between the keel and wave surface 

in degrees 

If the wave is assumed approximately sinusoidal the wave 
slope can be calculated as 

0 = 935. H(l/3) X / (16) 

where X = distance away from the slam center in feet 

Using equations (15) and (16) the most probable or the 
extreme pressure for design may be modified to allow for other 
than zero impact angle. 

The peak pressures do not occur at the same time along 
the length of the ship. From full scale measurements on the 
Wolverine State and model tests [14], it was determined that 
the pressure pulse traveled from forward to aft or vice versa 
at speeds from 260 to 520 fps for a 520 ft ship. Ochi assumes 
that the slowest traveling velocity is the worst slam condition 
and that the results can be Froude scaled for different length 
ships to get 

v^ = 260 /L/520 = 11.4 /L (17) 
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where 



= traveling velocity in feet per second 
This assumption was not felt valid since the traveling velocity 
should increase as the rate of entrance into the water in- 
creased — that is, as relative velocity increased so should the 
traveling velocity. It is agreed that if a slam could occur 
at the highest possible relative velocity and the slowest 
traveling velocity it would be the worst case, but it is felt 
that these two conditions cannot be met at the same time. 

The traveling velocity is calculated using the slam 
model in Figure 6. Using this model, the height above the 
water surface at the slam center, the most probable or extreme 
relative velocities of the slam center, and the length of the 
slam area are known or can be found. From this information 
the traveling velocity is given by 

''t * 

where 1^ = slam length from FP to aft extent of slam 

t. = time slam takes to travel from aft extent 
^ of slam to slam center 

and 



t^ = H/VEL 



(17b) 



where H = height of slam center above the wave surface 
at start of slam 

VEL = most probable or extreme relative velocity 
of the slam center 
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The last information required before integration can 
proceed is the spatial distribution of pressure along the 
girth. The most probable and the extreme pressure for design 
as calculated by equations (10) and (12) and modified with 
equation (15) for wave slope are only valid at the keel. 
Following again Ochi's work [14], [18], the maximum slam 

pressure can be assumed constant along the flat bottom of the 
hull. From the end of the flat the pressure drops off approxi- 
mately linearly to zero at the 1/10 design draft, T' . Figure 
8a shows this for a typical section. The pressure component 
of interest is the vertical component on the hull. The values 
in Figure 8a must be multiplied by COSB where 3 is defined as 
the angle between the normal and the vertical at the point in 
question. Figure 8b shows the new vertical component of the 
pressure distribution needed to find the slam force. 

All the information is now available to calculate the 
slam force and impulse. The procedure is as follows: 

1. Using equations (4) and (5) along with the r.m.s. 
relative motions and velocities from tables [15] , 
calculate the nu nber of slams at each station. 

The last station aft with N(t) - 1 is the longi- 
tudinal extent of slam. 

2. Calculate the most probable extreme pressure from 
equation (10) or the extreme pressure for design 
from equation (12) at each station forward of the 
last slam station as found above. 
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FIGURE 8a - GIRTH PRESSURE DISTRIBUTION 




FIGURE 8b - DISTRIBUTION OF VERTICAL PRESSURE 
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3. Assuming the carter of slam is midway between the 
FP and the last station where slamming occurs, 
calculate the slopes using equation (16) and 

the corrected pressures at each station using 
equation (15) . 

4. Calculate the vertical pressure components at 

each station for approximately 10 waterlines between 
the baseline and T' as in Figure 8. The spatial 
distribution of the vertical pressure is now known 
over the entire slam region. 

5. Assuming a triangular shape as in Figure 5, and 
a pulse duration of t^^ as predicted by either 
equation (13) or (14), add temporally at each 
location. The start of slam pressure rise (t^) 

at each spatial point is predicted by equation (17a) 
in the longitudinal direction and by the relative 
velocity in the girth direction. 

Due to the numerical integration procedure necessary, a 
computer program was written to perform the above calculations 
and numerical integrations. The program and the details of 
the above procedure are discussed in the next section. 

By evaluating the pressures at different time increments, 
the force time history can be determined for each section. 
Integrating these force histories over time will give the 
slam impulse at each station. 
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V. COMPUTER PROGRAM AND DISCUSSION OF RESULTS 



A computer program was written to implement the procedure 
discussed in Section IV. A listing of the program is included 
as Appendix B. There are many input and output options which 
are discussed in Appendix C which is designed to be a user's 
manual. A short example program illustrating the data deck 
and some of the possible outputs is included as Appendix D. 

In this section the results obtained from the program will be 
discussed using the Mariner and Fotini-L as example ships. 

Throughout this thesis stations are numbered from 0 at 
the forward perpendicular to 20 at the after perpendicular. 
Since slamming occurs in the forward portion of the ship where 
the lines change rapidly, half-stations were introduced num- 
bered from 0 to 40 similarly. The slam forces calculated are 
considered as centered at the half-stations extending equal 
distances forward and aft of these half-stations. 

Two computer programs were written and used. The assump- 
tions and calculations of Reference [14] were duplicated and 
used to compare with the results in this work. The only 
deviations were the use of equation (6a) for k > 0.19 and the 
use of seakeeping information obtained from the Reference [15] 
tables. The results for the Mariner were identical to those 
obtained by Ochi [14] . The significant aspect of this is 
that good results were obtained without recourse to seakeeping 
computer programs. 
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OPERATING CONDITIONS 



The Mariner was used as one of the example ships and the 
computer program was run inder various operating conditions. 
The four conditions selec :ed for discussion here are: 



1. 


MARINER 


FULLY 


LOADED 


7.4 


KNOTS 


2. 


MARINER 


FULLY 


LOADED 


15.5 


KNOTS 


3. 


MARINER 


LIGHT 


DRAFT 


15.5 


KNOTS 


4. 


MARINER 


LIGHT 


DRAFT 


7.4 


KNOTS 



all at a significant wave height of 25 feet. Table la lists 
the ship characteristics for the two loading conditions con- 
sidered. This ship and conditions were chosen since they 
were used in Reference [14]. 



MARINER 


FULLY LOADED 


LIGHT DRAFT 


BLOCK COEF. (C„) 


0 


. 613 


0.583 




L/B 


6 


.95 


6.95 




B/T 


2 


. 82 


3.80 




DRAFT FWD. 


27 


FEET 


16.4 FEET 


DRAFT AFT. 


27 


tt 


23.6 


tt 


H(l/3) 


25 


It 


25 


tt 


LBP 


528 


tt 


528 


tt 


BEAM (B) 


76 


tt 


76 


tt 


DESIGN DRAFT 


29. 


75 " 


29.75 


tt 



TABLE la - MARINER OPERAf ING CONDITIONS AND DIMENSIONS 
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A second ship, the F5TINI-L, was analyzed and is one of 
the three ships reported )n in Reference [4]. The principal 
characteristics and loadi ig conditions are given in Table lb. 



FOTINI-L 

BLOCK COEF. 
L/B 
B/T 

DRAFT FWD. 
DRAFT AFT. 
H(l/3) 

LBP 

BEAM 

DESIGN DRAFT 
SPEED (KNOTS) 



FULLY LOADED 

0. 84 
7.55 
2.79 

38.0 FEET 

38.0 " 

20 to 40 FEET 

800 FEET 

106 

38.0 " 

9.59, 14.38 



BALLAST 

0.65 

7.55 

3.55 

28.31 FEET 
31.46 " 

20 to 40 FEET 

800 FEET 

106 

38.0 " 

9.59, 14.38 



TABLE lb - FOTINI-L OPERi^TING CONDITIONS AND DIMENSIONS 



R.M.S. RELATIVE MOTIONS ;^ND VELOCITIES 

The first modification to the Ochi-Motter approach was 
to apply a broadness factor correction to the r.m.s. relative 
motion and velocity. If the ship responses are assumed to be 
narrow band processes, the r.m.s. values of the responses can 
be obtained directly fron the tables. These tabular results 
do not take into account the fact that seakeeping processes 
are not ideally narrow be nd [15] . More realistic values can 
be obtained if a broadness factor correction is applied: 
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(18) 



m 



(1 - e^/2) 



where = r.m.s. value of response from tables 

= corrected r.m.s. value 
e = broadness factor 



and 




th 

where m^^ = the moment of the spectrum. 

A value of e = 0.59 is used in the calculations [7], [15]. 

This correction factor reduces the relative motions and 
velocities and therefore reduces the slam forces predicted. 

RESULTS USING OCHI-MOTTEJ' METHOD 

The results for the four operating conditions using the 
Ochi-Motter method with the broadness factor correction are 
given in Tables 2 to 5 and Figures 9 to 12. These results are 
for each of the half-stations where slamming is predicted. 
Table 6 shows the total slam impulse and the amplitude and 
duration of an equivalent slam for the ship. Table 6 results 
are presented for ease of comparison of the two methods. 
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TABLE 2 ■ ^MARINER - OCHI METHOD 
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TABLE 5 - 4-MARINER - OCHI METHOD 



E30T 




I 




(SdlV) 33U0d 
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FIGURE 9 - 1-MARINF R - fiCMl MODEL 



E30T 




(Sdl^) 33dOJ 



53 



FIGURE IQ - 2-MARINER - OCHl MQJEL 




(SdIM) 33MDJ 



s-i 
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FIGURE 12 - 4-MARIN FR - OCHI MODEL 



MARINER 



TOTAL SINE FORCE APPROX, 



TOTAL IMPULSE 



CONDITION 


AMPLITUDE 

(KIPS) 


DURATION 

(SEC) 


(KIP-SEC) 


1 


985 


0.43 


269 


2 


4490 


0.62 


1766 


3 


4406 


0.68 


1901 


4 


1220 


0.48 


372 


TABLE 6 - 


TOTAL SLAM FORCE 


AND IMPULSE 


(OCHI METHOD) 



NEW METHOD RESULTS 

The changes made to the Ochi method discussed in the 
previous sections and incorporated into the second computer 
program are: 

1. an angular correction applied as per equation (15) , 

2. a different model of slamming as shown in Figure 6, 
and 

3. a traveling velocity calculated by equation (17a). 
The results for the same four operating conditions are given 
in Tables 7 to 10 and Figures 13 to 16. Table 11 gives the 
total values of slam force and impulse corresponding to 
Table 6 above. 
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TABLE 7 - 1-MARINER - NEW METHOD 
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TABLE 9 - 3-MARINER - NEW METHOD 
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TABLE 10 - 4-MARINER - NEW METHOD 
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FIGURE 13 - 1-MARINER - NEW MODEL 
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FIGURE 14- 2-MARINER - NEW MODEL 
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FIGURE 15 - 3-MARINER - NEW MODEL 
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FIGURE 16 - 4-MARINER - NEW MODEL 



mariner 



TOTAL SINE FORCE APPROX. 



TOTAL IMPULSE 



CONDITION 


AMPLI rUDE 
(KIPS) 


DURATION 

(SEC) 


(KIP-Sl 


1 


1632 


0.18 


187 


2 


7312 


0.21 


978 


3 


7105 


0.22 


995 


4 


2084 


0.17 


232 


TABLE 11 - 


TOTAL SLA 4 FORCE 


AND IMPULSE 


(NEW MODEL) 



COMPARISON OF RESULTS BETWEEN THE TWO METHODS 

In general the new model predicts lower slam impulses, 
shorter slam durations an i higher slam forces. In addition, 
the time histories are more triangular in the new model. 

The lower slam impul 3es are a result of the angular 
correction term only. If the correction were not used, both 
methods would yield the s me impulse values. 

The shorter total duration is a result of the slam model 
and the calculated traveling velocity. The model reduces the 
duration by approximately 1/2 since the pressure wave travels 
forward and aft toward the slam center. An additional decrease 
in slam duration occurs because the calculated traveling 
velocities are, in general, higher than those calculated by 
Ochi. For these four operating conditions the traveling 
velocities range from 731 fps for case 3 to 1016 fps for 
case 1 . 
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The total force on the ship is the spatial and temporal 
integration of the local pressures. In the new method the 
pressures are always less than or equal to the Ochi method 
but the new method yields larger forces. This is due to the 
integration over time. Referring to the plots of both methods. 
Figures 10 and 14, the reason for this can be seen. The half- 
station forces are more bunched together in the new method. 

The total force on the ship at any time is the sum of the 
forces on each half-station at the same time. Due to the 
greater overlap in the nev/ method, more half-stations are 
summed at the same time, resulting in higher forces. 

The importance of the various parameters of storm duration, 
threshold velocity, local slam duration, significant wave 
height, and most probable vs. extreme pressure for design 
were looked at to determine the sensitivity of the results 
to their variation. 

SENSITIVITY TO STORM DURATION 

To determine how important accurate prediction of the 
maximum storm duration was, the MARINER case 2 was run for 
different values of storn duration with all other input 
variables held constant. As seen in Figure 17 the total slam 
impulse is not sensitive to storm duration. A 50% increase 
in storm duration from 3C to 45 hours results only in a 4% 
increase in slam impulse. 
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FIGURE 17 - TOTAL SLAM IMPULSE vs. TIME OF STORM 
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SENSITIVITY TO THRESHOLD VELOCITY 



Since the threshold velocity at which slamming occurs 
has been defined differently by various authors, the sensi- 
tivity of slam impulse to changes in threshold velocity was 
determined. Values of threshold velocity between 6 and 20 
ft/sec were used for MARINER case 2 and the total slam impulse 
plotted against threshold velocity in Figure 18. From this 
figure it can be seen that the use of 12 ft/sec is reasonable 
and, at worst, conservative. 

SENSITIVITY TO LOCAL SLAM DURATION 

The total slam impulse on the ship is the sum of all the 
local impulses. Those parameters which affect local impulse 
also influence the total impulse the same amount. From 
Figure 5 the local slam impulse is given by: 

Ij = 1/2 p_^(o) (20) 

2 . 2 

where Ijj = local impuj se/in = Ib-sec/in 

From (20) it can be seen that a large variation in t^ 
results in a large change in Ip and subsequently in the total 
impulse. Therefore total slam impulse is very sensitive to the 
local slam duration. Since the determination of t^^ is in 
question, this is considered the area where further work is 
necessary (Appendix A) . 



68 




69 



I 



I 



I 



AFFECT OF SIGNIFICANT WAVE HEIGHT 



Another ship, the Fotini-L, was also analyzed. Several 
computer runs were made using both programs varying the signi- 
ficant wave height from 12 to 40 feet. Figure 19 is a plot of 
the results. The Ochi method gives results approximately 2 1/2 
times the new model results. As expected, the impulse in- 
creases with significant wave height and with velocity in the 
seaway. The rate of increase decreases at higher significant 
wave heights and is expected since the ship responses behave 
similarly at higher sea states. 

MOST PROBABLE vs. EXTREME PRESSURE FOR DESIGN 

The four Mariner cases were run using the new method for 
the most probable extreme pressure and the extreme pressure 
for design cases. Tables 11 and 12 summarize the results. 



MARINER 

CONDITION 


TOTAL SINE 
^sj^PLITUDE 
(KIPS) 


FORCE APPROX. 
DURATION 
(SEC) 


TOTAL IMPULSE 
(KIP-SEC) 


1 


733 


0.20 


91 


2 


3248 


0.25 


510 


3 


3 663 


0.25 


591 


4 


] 147 


0.19 


137 


TABLE 12 - 


TOT7L SLAM FORCE AND 


IMPULSE (MOST 


PROBABLE 




EXTREME PRESSURE USED) 
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mrm. 



The duration of the slam is greater in the most probable case 
since a lower relative velocity predicted results in a slower 
traveling velocity (n.b.. Table 11). The impulses are approxi- 
mately 1/2 those obtained from the extreme pressure for design 
calculations. Since the most probable extreme value is indi- 
cative of the largest value observed in a trial, these values 
should be used if any comparison is done with full scale data. 

SHAPE OF THE SLAM FORCING FUNCTION 

Knowledge of the shape of the slam forcing function is 
necessary to input the slam force into the S.H.V.R.S. program. 
The S.H.V.R.S. program will accept three input shapes for the 
forcing function; 

1. triangular 

2. half-sine wave 

3. damped half-sine wave. 

The forcing function can be applied at each of the half-stations 
with time delays corresponding to the traveling velocity or a 
single slam force can be applied at one location. Since much 
of the previous work [2], [3], [4] used a single sine pulse, 

this program develops a half-sine approxim.ation to compare with 
these results. 

The shape of the force at each half-station is different 
in the Ochi method and this method. In the Ochi method the 
assumed shape of the local slam (Figure 5) is not critical. 
Several pulses of different shapes were tried in the Ochi 
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program but the predicted half-station forcing functions were 
about the same. This is because the traveling velocity is 
slow enough that the slam starts at different times along the 
length of the half -station . In the new method the traveling 
velocity is, in general, high enough that the slam starts at 
the same time along the entire half-section. This results in 
a force shape the same as the assumed local slam pulse. The 
shape of the slam force at each half-station (Figures 9 to 12) 
in the Ochi method can be approximated as half-sine functions. 

The program calculates the amplitude and duration of these 
approximations. The shape of the force at each half-station 
(Figures 13 to 16) in the new method can be assumed triangular, 
of height equal to the maximum force, and of duration 

= 1.1 = 4.8 X lO"^ /l (21) 

where t 2 was determined by neglecting the "tails" of the 
force function and maintaining the same impulse. 

The total slam force shape can be obtained by adding the 
forces from each half-station at the same time. Figure 20 is 
a plot of the total slam force for the MARINER case 2, time of 
storm 30 hours, calculated by the new program. Figure 21 is 
the plot of the forcing function at each half —station . If the 
points which are less than 1% of the maximum value are neglected 
and a linear regression analysis done on the remaining points 
in Figure 20, the triangle formed is almost equilateral. There- 
fore this shape should be used for the total slam forcing function. 
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FIGURE 20 - TOTAL SLAM FORCE vs. TIME 
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FIGURE 21 - 2-MARINER - NEW MODEL - TSTORM r 30 hours 



midship slam bending moment and stress 



The FOTINI-L was considered at different significant 
wave heights, velocities, and loading conditions. The results 
for the total slam impulse in the ballast condition are shown 
in Figure 19. As the significant wave height and/or ship 
speed is increased, the slam center moves aft. Donovan [4] 
determined midship bending moment and stress for the FOTINI-L 
using a 100 ton-sec half-sine slam impulse, applied 36 feet 
aft of the forward perpendicular. This slam location is close 
to half-station 2 of this model. From the computer runs on the 
FOTINI-L, two operating conditions resulted in slam centers at 
half-station 2. These are given in Table 13. 



CONDITION 


LOADING 


SPEED 


H(l/3) 


TOTAL 


DURATION 


CONDITION 


(KNOTS) 


(FEET) 


IMPULSE 

(KIP-SEC) 


OF IMPULSE 
(SEC) 


1-FOTINI-L 


FULL 


14.38 


32 


1199 


0.18 


2-FOTINI-L 


BALLAST 


9.59 


20 


960 


0.18 


TABLE 13 - 


FOTINI-L 


- SLAMS 


PREDICTED 


AT HALF-STATION 2 



From the graphs of response in Donovan's work, these 
slam impulses give the moments and stresses in Table 14 for 
various stiffnesses. This is determined assuming a linear 
relationship between bending moment and impulse which holds 
if the slam durations are the same — that is, holding all else 
constant, doubling the slam impulse doubles the bending 
moment response. 
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CONDITION 


60% 


BENDING 


1-FOTINI-L 


343 


MOMENT 






(FT- TONS) 


2-FOTINI-L 


246 



STRESS 


1-FOTINI-L 


26.5 


(KPSI) 


2-FOTINI-L 


21.9 



STIFFNESS* 



80% 


100% 


120% 


140% 


364 


375 


391 


404 


268 


283 


296 


313 


22.1 


18.9 


16.5 


14.5 


19.2 


16.7 


14.5 


13.1 



TABLE 14 - FOTINI-L - BENDING MOMENT AND STRESS AMIDSHIPS DUE 
TO SLAMS CENTERED AT HALF-STATION 2 

*The stiffnesses are % of as built values. 



The values in Table 14 for the fully loaded case are 
approximately 5.4 times what was predicted in Reference [4] 
and 4.3 times the predicted results in the ballast condition. 
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VI. RECOMMENDATIONS FOR FUTURE WORK 



As seen by the values of stress obtained in Table 14, 
the slam forces and impulses result in large midship stresses 
on the FOTINI-L. The stresses could be evaluated only for 
slams at half-station 2 since dynamic response results were 
available only at this station. As the slam center moves aft, 
the total slam impulse and force increases. Higher slam 
stresses could result from slams at different locations. 
Therefore additional slam locations should be considered in 
the S.H.V.R.S. program in any study of the slam response. 

The shape of the total slam force is more triangular 
than sine (n.b.. Figure 20). The S.H.V.R.S. program should 
be run for several triangular pulses and the results compared 
to the half-sine results available. If they do not come close, 
a triangular pulse should be used. 

The validity of using a single total force approximation 
vice distributed forces along the bottom needs to be verified. 
The results of this work give the distributed forces as well 
as a half-sine approximation for the total slam force. Both 
should be run on the S.H.V.R.S. program to determine if the 
two agree. If they do not, either the distributed forces 
should be used or a new approximation developed which does 
agree . 

Although the slam impulses are less in this method than 
in the Ochi method, they still result in large slam stresses. 
There is still reason to believe, because of these large 
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values, that these impulse values are also overpredicted. As 
discussed earlier, and in Appendix A, equation (13) for the 
duration of the local slam pulse is in question. There pre- 
sently exists a large amount of literature on drop tests but 
few model and full scale results are available to determine 
durations. Correlation with full scale data needs to be 
accomplished to determine if equation (13) is valid or some- 
thing like equation (14) would be better. 

A series of model tests should be run to verify that the 
model of slamming used here is valid. This would include 
verifying or modifying the traveling velocity calculations 
and the correction equation for the angle between the keel 
and the wave surface (equation (15) ) . 
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VII. CONCLUSION 



The work presented here develops a slam force prediction 
method which is sensitive to seakeeping parameters based on 
a model of slamming different from other methods developed 
to date. The results obtained show slam impulses and time 
durations less than those in the Ochi-Motter method [14] , but 
impulses greater than the 100 ton-sec value used by Kline [3] . 
It is felt that values closer to those predicted here should 
be used in the overall A.B.S. study on hull stiffness. As a 
minimum, the results of this work have shown that the 100 ton- 
sec value used in the work to date is in question and should 
not be applied to all ships under every operating condition. 
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APPENDIX A - LOCAL SLAM PRESSURE, DURATION OF PULSE 



Ochi predicts that the slam pulse is triangular and the 
duration is 0.1 seconds for a 520 foot ship. From this he 
Froude sacles and concludes that equation (13) holds - i.e., 

= 4.4 X lo"^ yz 

From full scale measurements on the USCGC UNIMAC [19] , 

L = 300 feet, values of slam duration locally were of the 
order of 0.020 to 0.025 seconds. Using equation (13), this 
length ship results in a predicted duration of t^^ = 0.076 
seconds. This discrepancy was too great to ignore without 
further consideration. 

From drop tests done by Chuang [8] , [9] , and [20] , on 

flat bottom impact with zero deadrise angle, the following 
equation holds for deadrise angles less than 3® 



ti 4 






(22a) 



where 

1 = half-width of infinitely long flat plate 

o 

C . = speed of sound in air = 1125 ft/sec 

air 

For the UNIMAC and the durations of slams measured this 
would predict the half-width at frame 23 to be of the order 
of 6 feet. Lines were not available to determine how close 
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this agrees to the half-breadth of the 1/10 draft at this 
section, but it is a reasonable value for a ship with beam 
41 ft. 

In other drop test data [21] equation (22a) held quite 
well for many different sized test specimens if the beam was 
used as the scaling parameter. Without more information from 
full scale tests, equation (22a) cannot be verified and there- 
fore will not be used. It is felt that the local half-breadth 
should be the scaling parameter, not the length as used by 
Ochi . 

It is proposed here that 



which comes from equation (22a) would be a better estimate 
of slam duration with possible modification of the constant 
from full scale data. 

Equation (22a) holds for small deadrise angles less than 
3° where the compressibility of air influences the duration of 
pressure pulse. For larger angles the impact pressure time 
history becomes more steep fronted with about half the duration 
predicted by equations (22a) and (14) . Therefore above 3° 
equation (22b) holds fairly well. 



t 



1 



3.6 X lO"^ B* 



(14) 




(22b) 
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Although equations (22a) and (22b) are for deadrise angles, a 
similar result might be expected for angles between the Iceel 
and the wave surface. Therefore equation (22a) is probably 
valid in the flat region near the slam center while equation 
(22b) is probably valid outside this region. 

For the MARINER equation (14) yields values of t^^ = 0.022 
at half-station 4 to t^ = 0.115 at half-station 16. Therefore 
it results in much lower impulses than would result using 
equation (13). For the FOTINI-L larger values were found 
using equation (14) than from equation (13). At half-stations 
4 and 10, t^ = 0.130 and 0.189 respectively using (14). This 
would result in larger in pulses than one would get if (13) 
were used. If equation (22b) were applied in the appropriate 
regions, the above resulting impulses could be reduced. 

Use of equations (22a), (22b), and (14) are not justified 

at this time, but it is felt that full scale data should be 
analyzed to determine if equations of this form are valid. If 
they are, further reduction of calculated impulses by the 
program is possible in many ship cases. 
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APPENDIX B - COMPUTER PROGRAM LISTING 



This appendix contains a listing of the FORTRAN computer 
program. It was written for an Interdata Model 70 computer 
and requires approximately 60 K of memory. Two subroutines 
are called which are not listed here but are readily available 
to most users as part of the IBM Scientific Subroutine Package. 
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C PROGRAM SLAMCW WALKER JANUARY 1976 

C program main THIS PROGRAM CALCULATES THE SLAM FORCES AND 

C IMPULSES generated ON A SHIP IN AN IRREGULAR SEA 

DIMENSION XIN(5) •YIN (5) •XOUT(H) .YOUTdl) »XX ( 1 ?• 1 1 ) • YY ( 1 1 ) tT (20) • 

C VK(17) •YO(ll) •RM(17) •PV(17) fEPREScn) .PRnRPRd?) * 

C PRES <17* 11) •FO^CE( 17.20) .DVEL(17) .T IMF (17) 

DIMENSION A (17) .OUR (17) .SLAMS (17) .TM(20) •SHAPF(17) • 

CPRSLAM(17) .FMAX(17) •AIMPLS(17) .PVEL(17) ♦XG(17.11) .XSCL(A) 
data ESL*ESP.NIN.N0UT,PI/99..99..5. 11 •3.1<^1592654/ 

DATA N.M/R,S/ 

DATA shape /O. Q 9375. 0. 79 1 6667. 0. 6937S • 0. 6f 0.5 104167 •0.425.0* 3475. 
C0.266666*0.19375.0. 125.0. 1 6041 666. 0.00. -0.05625. -0.1 0R333. -0.15625 
Ci-0. 2. -0.239583/ 

1 READ (N.IOOO) 101. 102. 104. 108 

1000 FORMAT (7IS) 

IF (I02.LT.I) 00 TO 999 
C READ SHIP NAME AND loading CONDITION 

READ (N.lOOl) (T(I) .1=1.20) 

1001 format (20A4) 

C READ SHIP DATA CONSTANT IN PROBLEM 

READ (N.1002) RB.HH.ALBP*ALWL*STASP 

1002 FORMAT (5Fin.5) 

WRITE (M.1003) 

1003 format (d».35)(.*THIS PROGRAM CALCULATES THE SLAM FORCE OISTRIBUTl 
CON VS TIME*./.38X. *0N THE LOWER I/IO OF THF SHIP FROM STATION 0 FO 
CRWARD. •♦////) 

C WRITE SHIP NAME AND DATA 

WRITE (M.1004) (T (I) .1=1.20) .B0.HH. AL0P»ALWL.STASP 

1004 format (20X.20A4.//.30X. 'half BEAM = * .F10.5. lOX. ‘DESIGN DRAFT = • 
C»P10.5.//.30X.*LPP = ‘.F10.2.16X. ‘LENGTH LOAD WATER LINE s ‘.F10.2 
C.//f30 X. *DTSTanCE between stations = •.F6.2.40X.//.35X. • (all OISTA 
CNCES IN FEET) ».////) 

C CALCULATE THE K VALUES USING SUBROUTINE SPOO 
HsHH/10. 

S=HH/100. 

X0UT(1)=0.0 
00 10 1=2.11 
X0UTd)=X0UT(I-l) ♦ S 
10 YY(I)=X0UT(I) 

1*0 

C READ OFFSETS AT ANY 5 WaTERLINES 

READ (N.1005) ( X IN ( J) • J=1 *5) 

1005 format (5F10.5) 

20 READ (N.1005) (YIN( J) . J=1 *5) 

I*IM 

IF (YIN(5) .LE.0.5) GO TO 30 
C 

C CALCULATE OFFSETS IN 10 INCREMENTS FROM B.L. TO 1/10 DESIGN DRAFT 

c USING subroutine UGLYOK 

C 

CALL UGLYOK (NIN.NOUT . XI Ni YIN. XOUT ♦ YOUT .ESL^ESR) 

00 21 J=l.ll 

JJ=12-J 

Y0( J)=YOUT( JJ) 

21 XXd*J)=YOUT(J) 

B=YOUT(ll) 

IA0P=I01 

C CALCULATE THE FORM COEF. OF THIS STATION 

CALL SPOO (YO.VKK.B.H.S.BB.HH.IAOP.D 
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VK( I) =VKK 

IF (I04.lt. 1) GO TO 20 
0=FL0AT(I)/2. 

WRITE (M*1006) O, (YIN( J) t X I N ( J) t TOUT ( J) #XOUT( J) • Jal #5) * (TOUT (J) t 
C X0UT(J).J = 6*in 

1006 FORMAT ( M • *4QX.50 ( • ♦ • ) t / 1 SOX t • ♦ • ♦ 4«X • • ♦ • # /t SOX • • ♦ • • 1 1 X t *ST AT ION • • 
CFS.li* OFFSETS (FEET)*. 9X .**•./ .SOX* •*••♦ /.SOXt 50 (•••) ./• 
C50X. ,4PX* ••• */*SOX. ••• .9X. • INPUT • *9X. • • • # 7X . • C&LCUL ATEO • . 7X . 

C*/*50X* ••• *?3X. • .?4X. •«•./. SOX. SO (•••)./. SO X.*«*.lix.***,llx.*«* 
C.UX* ••• . 12X. ••• ./.SOX. *♦• .SX. *X* *SX. •♦•.SX. *Y* .SX. *♦• .SX. • X* .SX. • 
€•• *SX* *Y* .6X.*«*./.S0X.S0( •♦•) ./.SOX.***. 11X****.UX.***.UX.****1 
C2X*****/.S(S0X.***.F«.2.3X****.F8.2.3X. ***,FP.2.3X.***.F8.2.4X.*** 

C*/) *6(S0X, ***.1 IX •***•! IX.*** ♦F9.2.3X*****FR.2.4X****./).S0X.S0(** 
C*)*/**0*> 

GO TO 20 
30 NLINES=I-1 

READ (N.4003) IPLT. ( XSCL(K) .K=l .4) 

4003 format (IS.4F10.2) 

C CALCULATE AN AVERAGE FORM COEF. FOR FORWARD PORTION OF SHIP 

CALL KEFFEC (VK.AVGK) 

WRITE (M.lOll) AVGK . , 

loll FORMAT ( *0* .39X.S0 ( *** ) ./*40X. ••*.48X. *♦* */.40X. *♦* . 

c 6 x.*averagf hull shape Factor = *.f8.s.6x.***./. 

C40X. •♦* .4PX. •♦* • /.40X.S0 (•♦•)•//) 

C READ IN INTEGERS FOR CONTROL OF OUTPET OPTIONS 

100 READ (N.IOOO) 109. 103. 105. 106. 107. 1010 
IF (I09.LT.1) GO TO 1 
READ (N.1007) (T(I).I=1.20) 

1007 format (20A4) 

NSTA=NLINES 

WRITE (M.1020) (T(I) .1=1.20) 

1020 FORMAT (*1*.1X.20A4.//) 

READ (N.lOOfl) FROUOE. speed. SEA. HF.HA.C0.VAFT.VTHR 
1000 FORMAT (0F1O.S) 

READ (N.4000) TSTORM 
4000 format (FIO.S) 

HSIG=SEA*ALBP 

IF (TSTORM. GT. 0.01) GO TO SOI 
IF (HSIG.LE.3S.) GO TO SOO 
TST0RM=50.0-1.0^HSIG 
GO TO SOI 

500 TSTORM=78.0-1.8*HSIG 

IF (HSIG.GE.20.) GO TO SOI 
TST0RM=4S.-0.1S*HSIG 
501 CONTINUF 

IF (FROUDE.GT. 0.00001) GO TO 110 
FROUOE = SPFEn*1.6fifl9/SQRT(32.174*ALWL) 

110 SPEED = FR0UDE*S0RT(32.l74*ALWL)/1.6B89 
FSCALE=S0RT ( ALPP/520 . ) 

THETA s(HA-HF)/ALBP 

vafts=vaft*fscale 

vthrs=vthr*fscale 

ALB=ALBP/(?.*B8) 

ARTa4.*00/ (HA^HF) 

IF (lOS.EQ.l) GO TO 103 

READ (N.1008) (RM( I) .RV( I) . I=2.8.6) 

C CALCULATE THE RELATIVE MOTIONS AND VELOCITIES AT ALL HALF-STATIONS 

CALL RELHV (RH.RV) 

GO TO 104 
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103 NP=0 

122 IF (NR.GE.NSTA) GO TO 104 

RE AO (N, 1000) RMCNRol ) ,RV(NR^1) *RM(NR^2) *RV(NR^2) ,RM(NR*3) t 
CRV(NR^3) *RM(NR*4) ,RV(NR*4) 

NRsNR«4 
GO TO 122 

104 CONTI^4UE 

c calculate «;lam pressures at the keel* number of slams, extent of 

C SLAM, velocities OF SLAM. AND PROBABILITY OF SLAM 

CALL SLAMPP (RM. RV.HF, HA. ALWL. NSTA. ERRES.VK. SLAMS. PROBPR.'nVEL.VTHR 
CS.TST0RM.HSIG.PRSLAM.PVFL*STASP) 

NCENTR = NSTA/2 ♦! 

IF (VAFT.GT.0.01 ) GO TO lOS 
IF (NSTA.FO.O) GO TO lOS 

Hl»<l,-SIN( (PI/2.)-0.0714*FLOAT (NSTA) ) )*HSIG/2, 

IF (I06.EQ.O) GO TO 106 
TT*H1/PVEL (NCENTR) 

GO TO 107 

106 TT=Hl/OVEl. (NCENTR) 

107 VAFTS=FL0AT (NSTA) ^STASP/ (4.*TT) 

vaft=vafts/fscale 

105 IF (NSTA. EO.O.ANO. 103. EQ.l) READ (N.1069) JJJ 
1089 format (13) 

WRITE (M.miA) CB.FROUOEfSPEEOfAL6.ABT.SEA.HSIG.TSTORM.THFTA.HF.HA 
c.nlines.nsta 

1014 FORMAT ('O'.lRx.tTHE SHIP OPERATING CONDITIONS ARFM.//. 

C20X. 'BLOCK COEF. = *.F7.4./. 

C20X.*FROUOE NO. = • .F7.4 f 25X. • SHIP VELOCITY = ».F6.2.* KNOTS'./. 
C20X. 'LENGTH TO BFAM PATIO = '.F6.3./. 

C20X.'BEAM TO ORAFT RATIO = '.F6.3./. 

C20X, 'NON-OIMENSIONAL SEA STATF = ' .F7. 4. lOX. 'SIGNIFICANT WAVE HEIG 
ChT s '.F5.0.' FEET'./. 

C20X.'TImE of STORM = f.FS.O.' HOURS'./. 

C20X.'TRIM angle (RAOIANS) = '.F7.4./. 

C20X. 'ORAFT FORWARO = ' .F5.2.24X . 'DRAFT AFT s '.F5.2f/. 

C20X. 'NUMBER OF SHIPS LINES READ = '. 13. 13X. 'NUMBER OF HALF STATION 
CS CONSIOEREO = '.13./) 

WRITE (M. 1015) VAFT.VAFTS.VTHR.VTHRS 
1015 FORMAT (20X. 'TRAVELING VELOCITY = '.F9.1.' FT/SEC FOR 520 FT SHIP 

C ' f // f 

C20X. 'TRAVELING VELOCITY = '.F9.1.' FT/SEC FOR THIS SHIPt//. 

C20X. 'THRESHOLD VELOCITY = '.F6.2.' FT/SEC FOR 520 FT SHIP*.//. 

C20X. 'THRESHOLD VELOCITY = '.F6.1.' FT/SEC FOR THIS SHIP'.//. 

C20X.'THE PROBABILITY OF EXCEEDING THESE SLAM FORCES (ALPHA) = 0,01 
C'f//) 

IF (NSTA.EO.O) WRITE (M.3030) 

IF (NSTA.FQ.O) GO TO lOO 

3030 format (40X.*NO slam calculations performed SINCE THE'./.40X. 
C'NUMBER OF SLAMS IS LESS THAN 1 FOR ALL STATIONS'.//) 

IF (I06.EO.O) GO TO ISO 

C CALCULATE SPATIAL OISTRIBUTION OF PRESSURE ALONG GIRTH 

CALL PRESS (PROBPR.NSTA.PRES.XX.YY.SfXG) 

GO TO 63 

150 CALL PRESS (EPRES.NST A .pRES. XX. YY.S. XG) 

63 CONTINUE 

C INTEGRATE PRESSURE TO GF.T FORCES FOR EACH HALF-STATION 

CALL FORCES (FORCE.NSTA.TiME.OELTAT.PRESf XG.ALBP.DVEL.VAFTS.STASP. 
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C YYfPVEL# 106) 

CALCULATE half-sine APPROXIMATION WITH EQUAL IMPULSE AND DURATION 
AT EACH half-station 

CALL APPROX ( NSTA*F0RCF*FMAX*AIMPLS# AfOUR*OELTAT) 

IF (1010. EO. 2) GO TO 30?0 
IF (1010. FQ.l) WRITF(M.3012) 

TO 12 format (M».126(*»*)*/*lXf2(t«t,9x),*«*.14x**«»fl7Xf»»*f6(llXt*«M 
Cf/*1X*2(**».9X).»**,» PEL MOTION * EXT PRESSURE •♦•UX* 

C#* MAXIMUM ••#11X«** amplitude ••.11X*»* TIME **f/f 

C1X*«» half • FORM t,34( •«•),# NUMBER ♦ SLAM ♦ SLA 
CM ♦ OF SINE • slam • SLAM •••/t 

CIX*** STATION • COEF • REL VELOCITY • MOST PROP PRESS ♦ OF SLA 
CMS • FORCE • impulse • APPROX • DURATION • STARTS 

C/*lXf ••• *9X* »• (K) ♦ (NON-DIM) ♦ • *6X , • ( PS I ) • t 6X ♦ • * • # 1 1 X • 

C»* (KIPS) ♦ (KIP-SEC) • (KIPS) ••.2(« (SEC) *•)•/• 

ClXfl26(**M) 

IF (1010. EO.l) WRITE (M,3013) ( J* VK ( J) ,RM ( J) ,EPRES ( J) fSLAMS ( J) . 
CFMAX (J) •AIMPLS(J) *A(J) ,DUR(J) •TIME( J) tRV(J) •PRORPR( J) .J=1,NSTA) 

3013 FORMAT (17( IX* ,3x*I3»3X# • • • *2X *F6 . 3* • • • *4X #F6. 3*4X* # • t ,6X*F6. 1 . 
CSX* 

C*«» *2(1X*F9.0*1X*«*«) *?XiFR. 1*1X, •*t *2X*F8.0f 1X*««**2(3X*F6.3*2X* 
C«**)f/*1X*»*»*2(9X*»*») ,4X.F6.3*4X*****6X*F6.1*5X* •♦•*6(11X***«) * 
C/*lX*126(*^*) */) ) 

IF (lOlO.EO.O) GO TO 3015 
IF (I06.EQ.O) WRITE (M.3011) 

IF (106. EO.l) WRITE (M.TOU) 

3011 format (30X*»««* extreme PRESSURE USED IN FORCE CALCULATIONS ♦♦♦•* 
C//) 

3014 format (30X.*»** most PrOBaBLE PRESSURE USED IN FORCE CALCULATIONS 

C ♦•••*//) 

GO TO 3020 

3015 DO 3000 KMsI.NSTA 
TM(1 )=TIME(KM) 

DO 3006 KN=2*20 

3006 TM(KN)= TM(KN-l) ♦ OELTaT 
WRITE (M.3001) KM 
IF (IO7.EO.0) GO TO 3010 
IF (I06.E0.0) WRITE (M.3002) 

IF (106. EO.l) WRITE (M.3004) 

WRITE (M.3005) (TM( J) *F0RCE(KM* J) * J=1 *20) 

WRITE (M*3(^07) 

3010 CONTINUE 

WRITE (M.300fl) EPRES(KM) »PROBPR(KM) *VK(KM) 

WRITE (M. 3(^09) RM(KM) .RV(KM) *SLAMS(KM) *PRSLAM(KM) ,0VEL(KM) * 

C PVEL(KM) *FMAX(KM) .AlMPLS(KM) *A(KM) * OUR (KM) .TIME (KM) 

3009 FORMAT ( 30X . »NON-D IMENS IONAL RELATIVE MOTION = • *F7.4*/* 
C30X**NON-OIMENSIONAL RELATIVE VELOCITY = NF7.4*/. 

C30X* •NUMBER OF SLAMS DURING EXTREME STORM s •*F0.O*/# 

C30X# ‘PROBABILITY OF SLAM IMPACT = ‘.FS.S*/. 

C30X* •extreme RELATIVE VELOCITY s ‘.FB.U^ FT/SEC‘ */* 

C30X.‘MOST probable relative VELECITY = ♦•F8.1f‘ FT/SEC‘*//« 

C30X* ‘MAXIMUM VALUE OF SLAM FORCE = ‘*F10.1.‘ KlPS‘*/t 
C30X* • iMPlfLSF VALUE = ».FlO.O*‘ KIP-SEC‘./* 

C30X* ‘AMPLITUDE OF SINF APPROXIMATION = ‘*F10.1*‘ KIPS‘*/* 

C30X* ‘DURATION OF SLAM ImPaCT FORCE = ‘*F0.5*‘ SEC‘*/t 
C30X#‘TIME SLAM STARTS = ‘*FB.5*‘ SEC»*//> 

IF (I06.EO.O) WRITE (M*3011) 

IF (106. EO.l) WRITE (M.3014) 

3000 CONTINUE 
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3020 CONTINUE 

3001 rOQMAT (•l»*?9Xf40(»»M •/•30X,»***38X*»*»f/t30X**»*fIOX»*MALF STAT 
CION ••I3*IIX«««»*/,30X*»«'*3flX*»«» * /*30X»40 (•••)) 

3002 FOPMAT (30X, ••• • I 7X t • • • ♦ 3 X • ' EXTPEmR FQPCE » » 3X ♦ • • • f/*30Xf ••• *3X* »T 
CINE (SEC) • *3X* tSXf «FOP DESIGN* *SX, ••* ♦/♦30X* ••• ♦ I7Xi *7X»* (K 
CIPS) • f7X* ••• */O0X*40( ♦•' ) ./f TOX* ••• f I7Xf *20X* ••• ) 

3004 FOPMAT ( 30 X • • ♦ • • 1 7X • * • • t 3X ♦ • MOST PROBABLE •» 3X *•••♦/♦ 30X ••••* 3X T 

CIME (SEC) • f 3Xf •«' *5X» • FORCE • » SXt ♦/♦30X, ••• t 17Xt ••• tTX. • (K 
CIPS) ' *7Xf ' ** ♦/♦30X,40 ( ♦^>» ) f/tlOX* 1 1 7X» • • ♦ ♦ ?0 X t • • • ) 

3005 FOPMAT (30X,»<**,5X,F7.3,5Xf««*»6X*F8.U6Xf*«*f ) 

3007 FOPMAT (30X»40( *•• ) *//) 

3008 FOPMAT (30X* 'EXTREME PRESSURE = ••FlO.Of* PSI • f/f30Xf 'MOST PPO0ABL 
CE PRESSURE = •*F10.0»* PSI • t /* 30Xt 'FORM COEF* (K) =',F10*5t/) 

TIMPLS=0.0 
DO POO I=lfNSTA 

200 TIMPLS=TIMPLS ♦ AIMPLS(T) 

TSLAMTrTlME(NCENTP) ♦OUP(NCENTP) 

ttt=timpls 

TA c PI«TIMPLS/(2.»TSLAMT) 

LOCl=NCENTP 

WRITE (M,10I7) TAfTSLAMT 

lOlT format (20X.' average SLAM FORCE * ' fF0*Of ••SINE (PI«T/' •FG.At • ) 

C KIPS',//) 

WRITE (M,ini8) LOCI 

1010 FORMAT (25X, 'MIDPOINT LOCATION IS HALF-STATION* ♦ I4,///) 

TIMPLS = 0.0 
DO 202 I=I,NSTA 

202 TIMPLS=TIMPLS^ AIMPLS ( I ) •SHAPE ( I ) 

TB=pI»TIMPL5/(2.^TSLAMT) 

WRITE (M,10I9) TPfTSLAMT 

I0l9 format (20X, 'MODAL AVERAGE SLAM FORCE = » *F8,0t ••SINE (PI*T/' iF6.4 
CtM KIPS'//) 

WRITE (M,10I8) LOCI 
WRITE (M,I069) TTT 

1069 FORMAT ( '0 ', lOx, 'TOTAL SLAM IMPULSE = 'tFS.Of' KIP-SEC'*///) 

IF (IO3.EQ.0) GO TO 100 

C CALL PLOTTING ROUTINE FOP TIME HISTORY OF FORCES 

CALL PLOT (TIME, FORCE, 0ELTAT,NSTA,IPLT,XSCL) 

GO TO 100 
999 WRITE (MOOIO) 

1010 FORMAT CO',' END OF JOPM 
STOP 
END 
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SUBROUTINE UGLYOK ( NIN ,nOUT . XI N . Y IN « XOUT . YOUT .ESU . ESR ) 

THIS subroutine takes the OFFSETS AT EACH INOUT WATERLINE AND 
USES SPLINE-CUBIC INTERPOLATION TO OETEPMINE THE 10 EOUALl Y 
SPACED WATERLINE OFFSETS BETWEEN RL AND l/lO DESIGN DRAFT 
REAL'B H»D.A«S.AE«BE.CE.DE,HALFf TwO«SIXfRAD,SLP.DSIN»OCOS 
dimension <IN(S).Y1N(S) .XOUT (in .YOUTd n .H(S) »D(S) , A(25) .S<5) . 
C AE(5) .8E(5) .CE(5) .DECS) 

DATA HALF/0. 5D00/. TWO/2. ODOO/, SIX/6. 0000/. RAO/1.74532R25O-0?/ 

NMl=NIN-l 

NM2=NHl-I 

NE0=NM2 

DO 1 N=1.NM1 

H(N)=XIN(N»n-XlN(N) 

1 0(N) = ( YIN(N»n-YIN(N) ) /M(N) 

IF (ESL.GT.RO.O) go TO 2 
NEO=NEO« 1 

2 IF (ESR.GE.qo.O) go to 3 
NE0=NE0»1 

3 NS0=NE0**2 
DO 4 N=1 .NSO 

4 A (N)sO.O 
J = 1 

L*1 

IF (ESL.GE.RO.O) GO TO 6 
A(l)=TWO*H(l) 

A(2)=H(1» 

SLP=ESL*RAD 

S«1 ) = (0( n-(OSIN(SLP>/DCOS(SLPI ) ) *SIX 

J=J«1 

m=l«nfo 

A(M)=H(1) 

L=H.I 

6 DO 5 N = 1.N<*2 

IF (N.GT.l) A(L-1)=H(N) 

A (L) =TWO* (H(N) .H(N.] ) ) 

IF (N.LT.NM?) A(L*1)=H(N»1) 

S(J)=(0(N»1)-0(N) )*SIX 
J=J«1 

5 L=L»NEO»l 

IF (ESR.GE.90.0t GO TO 7 
A(L-l) =H(NM1 ) 

A (L) =-TW0«H(NM1I 
L=L-NEO 
A(L)=-H(NM1 ) 

SLP=ESR*RAD 

S(J)=(D(NM1)« (OSIN(SLP) /OCOS(SLP) ) )*SIX 

7 CALL SIMO (A.S.NEO.kERROR) 

IF (ESL.LT.90.) GO TO B 
DO 9 N=1 .N*^2 

M=NM2-N*2 
9 S(M)=S(M-n 
S(1)=0.0 

8 IF (ESR. GT. 90.0) S(NIN)=0.0 
DO 10 N=l.NMl 

AE(N)=(S(N»1)-S(N))/(SIX*H(N)) 

BE (N) =HALF«S (N) 

CE (N) =0 (N) -H(N) * (TWO»S(N) ‘SIN*! > > /SIX 
10 0E(N)=Y1N(N) 

J=1 
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DO 11 N=l,NOUT 

ir (XOUT (N) ,6E.XIN(2) ) GO TO 14 
J=1 

60 TO 15 

14 ir (XOUTIN) .LT.XININHI ) ) GO TO 16 
JsNHl 

60 TO 15 

16 IF (XOUT (N) .GF.XINC J*1 ) ) GO TO 12 

15 Hl=XOUT (N) -XIN ( J) 

H2=Hl**2 

h3=h1*H2 

YOUT(N)=AE( J)*h3»BF( J)*h2«CE(J)*H1»D£(J) 
GO TO 11 
12 J=J*1 

60 TO 16 
11 COMTINUE 
RETURN 
END 
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SUBROUTINE SPOO < X , VKK . R * H, S, BB . HH , I AOP. I ST A ) 

SUBROUTINE SPOO ( X , VKK ,B»h,S ,BB .HH, I AOP, ISTA ) 

THIS program CALCULATES THE EORm COE^. K DEVELOPED BY OCHI TO RE 
USED IN predicting THE SLAM CHARACTERISTICS OE A SHIP ERQm THE 
SHIPS LINES. REEERENCE SNAME TRANSACTIONS 1R71 

DIMENSION X(ll ) ,XM(U ) .VH( in .ZH(11 ) .ROOT(S) .W(A) , 

2 xN(i 1 ) ,7 N(i I ) .XI ( in .x2(i 1) »2i (1 n .22(1 n 

2 ,WW(6) .ROOTR(S) .ROOTl (5) .COL(ft) 

DATA PI/3. lAl592f>/.EPS/0.0005/.N/R/,H/5/ 

IE (lAOP.EO.n WRITE(M.IOOl) B.H.B9.HH 
IE (lAOP.EO.l) WRITE(M.lon) (X(I).I = l.in 
IE (lAOP.EO.n WPITE (M.1002) ISTA 
EORMAT (2PX.* male-station number «.I3./) 

CALCULATE AREA AND MOMENT OE INERTIA BY SIMPSON'S RULE. 



A=0. 

DO 10 1=2. R. 2 
10 A=A. A.*x ( 1 ) »2.*X ( I ♦ n 

A = 2.»S/3.*(A»X(l)*X(in«‘*.*X(lO)) 

T1=0.0 

Al=0.0 

DO 20 1=2. R. 2 
AI = AI* 1 . 

T1=TI.A.*X(I)*(AI*S)**2 ♦ 2.*X(I.1)*( (AI*1.)*S)»»2 
20 AI=AI»1. 

TI= 2.»S/3.*(TT.A.*X(10)*(9.»S)**2*X(ln*(10.*S)**2) 
determine a1 .A3. as. and U. 

0L=H/8 

IE (ABS(Dl-1.) .L£. EPS) 0L=0L‘EPS 

A=A*2./PI 

C2=T1*B./PI 

TK1=B**2» (1 ..DL*DL**2) 

TK2=B*( 1 .-OL) 

TK3=-3.*8*(1..DL) 

TK4= ( A.TKl) /TK2 
TK5=TK3/TK? 

TK6=2./TK2 

01=DL**2.'».»0L-7. 

02=-0L**3.DL»*2 -3.*DL»3. 

D3=-DL‘5. 

0A=0L**3»0L**2. 3.*OL«3. 

D5=-0L»*2-9.»0L-12. 

D6=3.»DL*7. 

07 = -l .75*DL**^-0.5*OL»*3-OL**2-1 .SROLi-l .25 
W(n=-A."TKR 

W(2) =2. *R*TK6*D3-A.»TK5-6. 

W(3)=B*(2.*(DG.T''5»03).B“TK6»Dn-A.*TK4 

W(4)=B*(2.*TK4*03.B*(nS.TK5*Dl »0 .5*B*TK6*02) ) 
W(5)=B**2»(T><4»01 »0 .S»“* (3.*D4*TK5*02) ) 

V(6) =0 .5»B®*3* (T<4»D2.B*D7) -C2 
IE (lAOP.EO.n WRITE (m.3000) (W(I).I=1«G) 

3000 EORMAT (' w= • .ME15.5.3X) »/) 

MM=5 

DO 30 1=1. R 
J=7-I 

WW(1)=W( J) 
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CONTINUE 

CALL BAIRS (WWtCOL«MM,ROOTR,ROOTI.KK) 

ERROR=0. 000001 
NR = 0 

DO 40 1=1 ,5 
ROOT(I)=0.0 

ir (ABS(ROOTI ( I ) ) .GT.ERROR) 60 TO 40 
IF (ABS(R00TR( I ) ) .LT.FRROR) 60 TO 40 
MR=NR» 1 

ROOT (NR) =ROOTR ( I ) 

40 CONTINUE 

IF (NR.FO.l) 60 TO 3R 

DO 38 J=1 .3 

KB=J»I 

DO 37 I=KB.NR 

IF (ROOT ( J) .GT.ROOT( I) ) GO TO 37 
R=ROOT ( J) 

ROOT ( J) =ROOT ( I) 

ROOT ( I) =R 

37 CONTINUE 

38 CONTINUE 

IF (NR.F0.4) 60 TO 39 
IF (ROOT(4) .6T.R00T(5) ) GO TO 39 
R=ROOT (4) 

ROOT (4) =ROOT (5) 

ROOT (S)=R 

39 CONTINUE 

IF (UOP.EO.l) WRITE (M, 3001 ) ( ROOT ( I ) . I = 1 i 5) . NR, Kw 

3001 format (• ROOT = • »5(ri0.5,3X) ,/, • NR = •,I5,» KK = »,I5,/) 

IF (KK.GE.D GO TO 99P 

DO ?5 1=1 .N9 
U=ROOT ( I ) 

ALFA=2.»U/9 
A5=TK4/U*TKS,TK6*U 
Al= ( 1 .-OL) /ALFA-AS 
A3= ( 1 . ‘OD /ALFA-1 . 

IF (lAOP.FQ.l) WRITE (“,3002) Al,A3,A5 

3002 format (• A1 = • ,F10.S,3X,*A3 = • »F 1 0 .5, 3X, • AS = '.FIO.S,/) 

IF (A3 .LT. 0.) GO TO 2S 

IF (ABS(Al-FPS) ,LE. 1. .ANO. ABS(A3-EPS) .LE. I. .AND. 

2 ABS(AS-EPS) .LE. 1.) GO TO 26 

25 CONTINUE 
GO TO 998 

COMPUTE K-VALUE 

26 TK=-3. 59893 ♦ 2.41988*AI - 0.872855*A3 ♦ 9.62395*A5 
VK=EXP(TK) 

VKK=VK 

IF( lAOP.GT.O) WRITF (M.2000) A 1 , A3 , A5 »U , VK 
IF (lAOP.GT.O) WRITE (M.2004) 

A=A»PI/4. 

IF (lAOP.EO.l) WRITE (M.3005) A 
3005 FORMAT (lx.» half-section AREA OF l/IO DRAFT = '.FIO.S,/) 

IF (VK.LF.0.19) go to 109 
VKK= 0.027»B»B/A 

IF (lAOP.EO.l) WRITE (M.3004) VKK 
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3004 FORMAT (lX.« THE REGRE«?';I0N ANALYSIS NOT USED K = • 

Ctri0.6OX,i where K = 0.027 *B«0/A»,/) 

109 CONTINUE 
GO TO 100 
99a WRITE (M.2no2) 

100 CONTINUE 
RETURN 

1001 format ( IHI ,5RX, 12h»«* SPOO *** //28x , 72HTHIS POQGRAM COMPUTES T 

2HE conformal mapping COFFFICIENTS ANO The ><-VALUE/?8X. 

2 63HNEEOED TO COMPUTE SLAMmING PRESSJPE BY THE RFGOEs«;ioN EOUATION 
2. //2ax. 14HHALF BREADTH = .F10.5«5H FE«'T.3X, 

2 25HHEI6HT ABOVE The BOTTOM = ,F10.5»5h FEET./2BX. IAHmalF Bs’AM 
2= .F10.5«5H FEET.21 X.TmoRafT = .F10.5.5H FSFTl 

lOU format (/ 2BX.31HTHE VALUES X .X X ARF . F 1 0.S/40X . 1 H0» 

2 3X.IH1.7X.2HI0.SX, 9(F10.5/5RX» .FIO.S.SH FEET) 

2000 format (///2ax,74HTHE COEF"^ 1C IFNTS OF THE TRANSFORMATION OF THE C 
2IRCLE IN THE ZFTA-PLANE ARE//?ax,SHA = .ElS.a.xx.SHA = .E15.8, 

2 7X.5HA = .F1S.8/29X . IHl ,2SX. 1H3.26X, IHS///28X. IRHThf SCALE FACTO 
2R IS.41X. 14HTHC W-VALUF IS//2ex,4HU = , E 1 5 . 8 , 38X , 4HK = .E15.8) 

2002 format (///28x,GfHA SOLUTION FOR U HAS NOT BEEN FOUND. olfaSE CH 
2ECK THE INPUT DATA.//) 

2004 format (/36X .S7HWHERF < = Exo (- 3.599 ♦ 2.420 a - 0.873 A *9 
2.624 A ) /69X, IHI , lOX, 1H3. lOX. IMS) 

END 
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SUBROUTINE SLAMPR ( RM ,R V »HF • HA. ftLWL » I ST A » EPRES « VK , SL AMS tPRORPR * 

C OVEL . VTHRS.TSTORM.HSIG.PRSLAM.PvFL .ST flSP) 

THIS SUBROUTINE CALCULATES THE EXTRE'^E PRESSURE AT EACH STATION 
WHERE The pROPABILITY OF SLAM IMPACT IS GREATER Than 1 
DIMENSION RM( 17) ,Rv ( 17) .£PRES( 17) .VK ( 17) .PROBPR( 17) .DVELI 17) « 
CPRSLAM(17) .SLAMS (17) .PVCL(17) .X(17) 

DATA ALPHA.G.pi/0.01 .32. 17A.3. 1A1592S/ 

DATA M.N/S.B/ 

RSTARS=VThpS««2 
AL = 0.<?83*ALWL 
T=TSTORM 
NSTA=1STA 
DO 10 I=1.NSTa 
H= ( (ha-hf ) /AO . ) “Float ( i ) *hf 
ORM=2.*< (al*RM( I ) “O.BD) «*2) 

DRV=a.*G*AL« (RV ( I ) *0,R3) «*2 
PRSL am ( I ) = fXP (-H“H/0RM-RST ARS/DRV) 

SLAMS(I)=(3<^00.»T/(2.»PI) )»S0RT(DRV/3RM)*PRSLAM(1) 

IF (SLAMS(l).GE.l.) GO TO 11 
1STA=I-1 
N5TA=ISTA . 1 

WRITE (M.loo?) NSTA.TSTORM.SlamS(I) .HSIG 
1002 format (‘o*.* the numbpp of slams at half station ‘.id.* during a 

CSTORM OF DURATION «.F3.0.* HOURS = ♦.F7.4./.» FOR A SIGNIFICANT W 
CAVE height of ‘.FS.O.* FEET’.//) 

GO TO 100 
11 EX=1./SLAMS(1) 

EPRFS (1 ) =V< ( I ) « (RSTaRS-ORV*ALOG( 1 .-( 1 .-ALPHA) **FX) ) 

PRORPR ( T ) =V«' ( I ) * (RSTARS ♦ DRV*AL0G(SLAMS ( I ) ) ) 

DVEL ( I ) =S0RT (EPRES ( I ) /VX ( I ) ) 

PVEL ( I ) =SORT (PROOPR ( I ) /M< ( 1 ) ) 

10 CONTINUE 
1STa=NSTA 

WRITE (M.ionl) ISTA.SLAMS(ISTA) .TSTORM.HSIG 
1001 FORMAT (•!•.♦ additional STATIONS SHOULD BE CONSIDERED SINCE AT HA 
CLE station *.I3./.’ the NUMBER OF SlAMS = ’.FS.!.* FOP A STORM OF 
CDURATION ».F5.0.‘ HOURS’*/** THE SIGNIFICANT WAVE HEIGHT = *.F5.0. 
C’ FEET’.//) 

return 

loo CONTINUE 

NCENTRriST A/2*l 
HSTaSP=STASP/2. 

XCENTR=FLOaT (NCENTR) *HSTASP-HSTASP/2. 

DO 1 J=l. NCENTR 

X( J) =XCENTR- (Float (J) *HSTASP-HSTASP/2. ) 

1 CONTINUE 
J=NCENTR* 1 

DO 2 I=J.NSTA 
X(I)=X(I-1) *HSTASP 

2 CONTINUE 

DO 5 1=1 • ISTA 
A=-I12.*HSIG*X(I)/(AL**2> 

EPRES(D=EppES(I)*(10.“*a) 

PROBPR ( I ) =PRORPR (I)*(10.**A) 

5 CONTINUE 

RETURN 
END 
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subroutine: relmv (rm,rv) 

THIS subroutine takes THE INPUTEO RELATIVE MOTIONS AND VELOCITIES 
AT stations 1 K A AND ASSUMES LINEAR VARIATION TO OETERMINE VALUES 
AT THE 17 hale-stations 
DIMENSION RM(17) ,RV(17) 

SL0PE=(RM(P)-RM(?) )/6. 

B=RM(?)-SL0PE*P. 

DO 10 1=1.17 
RM( I ) =SL0PE*EL0AT ( I ) *8 
10 CONTINUE 

SLOPE = (RV(8)-RV(?) )/A. 

B=RV(?)-SLnPE*P, 

DO ?0 1=1.17 

RV (I ) =SlOPE«ELOAT ( I ) ♦ R 

20 CONTINUE 
return 

END 
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subroutine forces (FORCEfNSTAfTIME*OELTATtPRES»XX*ALBP»OVELfVAFTS* 
C STASP*TY*PVEL* lOiS) 

THIS SUBROUTINE DOES THE SPATIAL AND TEMPORAL INTEGRATIONS AT EACH 
HALF-STATION OF THE PRESSURES TO FIND THE SLAM FORCE TIME HISTORY 
DIMENSION FORCE (17,20) ,PRES ( 17, 1 1) ,XX ( 17t 1 1 ) tTIME ( 17) ,0VELU7) , 

C YY(17) ,PVELC17) ,r(lO) , FOR (17, 20) 

TSLAMrO. 1«SQRT ( ALBP/520,) 

0ELTAT=TSLAM/10, 

N0ELT=InT(STASp/(2.*VAFTS*0ELTAT) ) 

IF (NOELT.lt. 1) NOELTsl 
T»STASP/(VAFTS^2.) 

NCENTR=NSTA/2^1 
TIME(1)=0.0 
00 4 I=2*NCENTR 
TIMECI)=TIME(I-1)^T 

4 CONTINUE 
I^NCENTR*! 

00 s j=i*nsta 
TIME(J)=TIME(J-1)-T 

5 CONTINUE 

00 20 I=1,NSTA 

IF (I06.E0.1) GO TO 50 

OELY=OELTaT*OVEL(I) 

GO TO 51 

50 0ELY*0ELTAT^PVEL ( I ) 

51 FORCE (I*1)=0.0 
DO 21 J=2,5 

21 F0RCE(I,J)=(FL0AT(J-1)/5.)«(PPES(I*1)*XXCI,1) ♦ PRES C T ,2) • (XX C I *2) 
C-XX(I,1) ) ) 

00 22 J=6,10 

22 F0RCE(I,J)=(FL0AT(11-J)/5.)*(PRES(I,1H^XX(I,1) ♦ PRES ( I t 2) • ( XX ( I t 2 
C)-XX(1,1) ) ) 

00 28 J=ll,20 
28 FORCE ( I, J) =0.0 
00 23 J=3,ll 
F(1)=0.0 

• 0ELX*(XX(I,J)-XX(I,J-1) )*PRES(I*J) 

DO 24 JJ=2,5 

24 FCJJ)=(FL0AT(JJ-l)/5.)«nELX . , 

00 25 JJ=6,10 ’ • . 

25 F(JJ) = (FLOAT (11-JJ)/5.)*0ELX . ; 

00 30 L= 1*20 

M»INT (YY ( J-1 ) / (0VEL( I ) •float (L) •OELTAT) ) 

INCsL 

IF (M.LE*1) go to 31 

30 CONTINUE 
GO TO 20 

31 M»INC 
MM=M^9 

IF (MM.GT.20) GO TO 100 
00 26 JJ=1,10 
K*M^JJ-1 

26 F0RCE(I*K)=F0RCE(I*K) ♦ F(JJ) 

GO TO 23 

100 IF (M.GE.20) GO TO 20 
MM=20 
K*0 

00 27 JJ=M,MM 
K*K^l 
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27 FOPCE(I.JJ)=FORCE(If 
23 CONTINUE 

00 29 J=U20 

FORCE(I.J)=(STASP«rORCE(It J)/lOOO.)*U4./FLOAT(NOELT) 
29 FOR(I*J)=FORCE(If J) 

20 CONTINUE 

IF (NDELT.LE.n RETURN 
00 60 I=lfNSTA 
NEXT=1 
LAST*0 

00 70 M=2*N0ELT 
NEXT*NEXT*1 
LAST=LAST^1 
00 80 J=NEXT*20 
JJ=J-LAST 

FORCE ( I f J ) rFORCE ( I • J ) ♦ FOR ( I • J J ) 

80 CONTINUE 
70 CONTINUE 
60 CONTINUE 
RETURN 
END 
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SUBROUTINE PRESS (EPRES.NSTA.PRES.XX.YY.S.XG) 



THIS subroutine determines the spatial distribution OP the vertical 

COMPONENT OF the PRESSURE AS A FUNCTION OF DISTANCE FROM THE CENTERLINE 

DIMENSION FPRESd?) iPRFSdT.lI) . XX ( I 7. 1 1 ) . YY ( 1 1 ) tCTHETA(lO) .X(1I> 

CtXG( 17. I 1 ) 

DO 10 IsI.n'JTA 
PRES(1.1)=EPRES(I) 

X(l)=XX(I.l ) 

DO 12 J=1.10 

DELTAX = XX ( T . J‘ 1 )-XX ( I . J) 

HYP=S0RT (S**2*nELTAX*»?) 

X(J»1)=X(J)*HYP 
12 CTHf TA ( J) =0ELTAX/HYP 
DO 13 J=1 . 1 1 
13 XG(1.J»=X(J) 

OELTAPrEPPES(I)/10. 

DO J=2,ll 
A=FL0AT( J-1 ) 

lA PREsn . J) = (EPRES(I)-A*0E1-TAP)*CTHETA<J-1) 

10 CONTINUE 
RETURN 
END 
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SUBROUTINE KEFFEC (VK.AVGK) 

THIS subroutine calculates THE AVERAGE OF THE FIRST II HALF- 
STATIONS called the average form COE*"* 

DIMENSION VK(17) 

AVGKrO.O 

DO 20 J=l.ll 

AVGK=AVGK*VK(J) 

AVC.K = AVGK/1 1 . 

RETURN 

END 



subroutine approx (NSTA,FORCE.FMax.aIM?lS.A.OUR.DFLTAT) 

THIS subroutine CALCULATES FOR EACH HALF-STATION THE IMPULSE BY 
INTEGRATING THE FORCE TIME HISTORY AND COMPUTES A HALF-SINE 
approximation WITH TH<^ SAME IMPULSE ANO DURATION 
DIMENSION FORCE (17.?0).FMAX(17) .AIMPLSC 17) . A( 17) «nuR(17) .v(20) 
DATA PI/3. 1 A15R2<S54/ 

DO 10 I=1.NSTA 
DO 11 J=10.?0 

IF (FORCE ( I . J) .GE.O.OS) GO TO 11 
DUR( I) rFLOAT ( J-1 ) *DELTAT 
NPTS=J 
GO TO 12 
11 CONTINUE 

OU«( 1 ) =20, *f)ELTAT 

12 DO 13 J= 1.20 

Y( J) =FORCE ( I . J) 

13 CONTINUE 

CALL MAXFOR (Y.AMAX) 

FMAX ( I ) =AMAX 

CALL OSF (0ELTAT,Y,Y.20) 

AIMPLS(I)=X(2P) 

Ad) = PI*AIMPI.S (I ) /(2.»0UR( I ) ) 

10 CONTINUE 
RETURN 
END 
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SUBROUTINe MAXFOB (Y.AMAX) 

THIS subroutine determines the maximum force from the force time 
history at each half-station using a Parabolic fit through the 

HIGHEST 3 POINTS. 

DIMENSION v(20) 

ELAST = Y(?) 

DO 10 1=3, fl 
N=I 

IF (Y( n .LT.FLAST) GO TO 11 
FLAST = Y(I) 

10 CONTINUE 

AMAX = I.OEIO 
11 AA=Y(N-2) 

BB = Y (N-n 
CC=Y(N) 

B=0.5« <A.»PR-3.*AA-CC) 

C = AA 

A=BB-AA-B 

X=-P/(?.*A) 

AMAX=A*X»*2 ♦ B*X ♦ C 

RETURN 

END 



subroutine plot (TIME. forces. DELTAT, NSTA, IolT.XSCI ) 

THIS SUBROUTINE USES THE M.I.T. J.C.F. PICTR PlOTTING ROUTINE 
- IF THE PROGRAM Is TO BE USED ON ANOTHER MACHINE, THIS IS TwF 
ONLY subroutine that REO'JIRES CHANGING ALONG WITH THE RELATED 
CARDS IN MAIN 

INTEGEP«a IXPLOT ( 17 ) .XLaB(AO) 

DIMENSION PI TR( 34,20) .TIME (17) .FORCES ( 1 7 , 20 ) ,XSCL(4) 

DATA IXPLOT/IR, 1R,20,21 , 22 , 23, 24 , 25. 2S. 27 , 28 ,29, 3P . 31 , 32 . 33. 34/ 

DO 10 I=1.NSTA 

T=TIME(I) 

DO 20 J=1.20 

PLTR ( I . J) =EORCES( I , J) 

I 1 = 1*17 
PLTR( 1 1 , J) =T 
20 T=T*0ELTAT 
10 CONTINUE 

CALL PICTR (PLTR,34,XLAB,XSCL,NSTA.20,-1,-1,3,IPLT,0,1,IXPLOT) 

PAUSE 

RETURN 

END 
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noon 



subroutine BAIRS(ACOF.a,NORO,ROOTR,ROOTI,IER) 
real ACOF ( 1 ) «A ( 1 ) ,R00TR ( 1 ) »R00TI < 1 ) 
data T0L2/1 .OE-3/.TOLX/S.OE-6/ 

DATA TOL*ITMIN, I TMAX/ 1 , OE-S. 1 0 ♦ 1 00/ 

DATA TOL3/1 .OE-10/ 

N=NORD 

IF (N.LT.0)N=-N 
0=AC0F (N» 1 ) 

OOM=-100.0 
IR = 0 

C BASIC ERROR CHECKING 

IF(N)1,1,2 

1 IR=1 

GO TO 200 

2 IF (N-<59)4,4,3 

I IR=2 

GO TO 200 

4 IF(0)6.5*6 

5 IR=4 

GO TO 200 

6 00 1 0 I = 1 , N 

10 A(I)=AC0F(I)/0 

A(N«1)=1.0 

C- REMOVE trivial zeroes FIRST 

J-2 

DO 14 1=1, N 
IF(A(I))12,n.l2 

II R00TR(J-1 ) =0.0 
ROOT I ( J- 1 ) =0 , 0 

14 J=J*1 

c FIND A REAL ROOT IF REDUCED POLYNOM lAL IS ODD 

12 LEFT=N-J»1 

IF (LEFT) 100,90,13 

13 NN=N»1 
IJ=J-1 

c estimate the modulus of The largest root 

DO 151 I=IJ,NN 
IF(A(I) ,NF,0,0) 

lOOMrAMAXl ( Al OG ( ASS ( A ( I ) /A ( IJ) ) ) /FLOA T ( I - IJ* 1 ) ,OOM) 

151 CONTINUE 

DOM=EXP(OOM| 

OJ=l ,0/OOM/DOM 
PJ=-2.0/OOM 

C SCALE the POLYNOMIAL SO LARGEST ROOT IS NEAR 1,0 

DO 9 I=IJ,N 
DO 9 K=IJ,I 
9 A(K)=A(K) /DOM 

IF(LEFT,NE.LEFT/2*2)G0 TO 15 
300 X0=0,0 

P0=A( J-1 ) 

C-- determine if real root is POSITIVE OR NEGATIVE 

Xl = l ,0 

IF(PO,GE,0.O)Xl=-l,O 
IP=1 

-EXPAND RANGE UNTIL ZERO HAS BEEN BRACKETED 

-THIS IS AN internal SUBROUTINE TO EVALUATE A REDUCED POLYNOMIAL 
400 ICNT=N*2 
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P1=0.0 

401 P1=P1»X1«A(ICNT-1) 
ICNT=ICNT-1 

ir (ICNT-J)40?.401 .401 

402 GO TO (320,330) , IP 
C 



320 


IF(P0*P1)325,360,321 




321 


xo = xl 

X1=10.0«X1 
GO TO 400 




c 


--ANSWER LIES BETWEEN XO AND Xl 




325 


IP=2 

X2=X1 




326 


X1=(X0.X2) »0.S 






IE(A8S(X2-XO)-ABS(X1*TOLX) ) 360, 


360,400 


330 


IF (P0*P1 ) 331 ,360,332 




331 


X2 = X1 
GO TO 326 




332 


X0 = X1 
GO TO 326 




c 


— NOW DIVIDE OUT THE ZERO 




360 


ROOTR(J-l)=xl*nOH 
ROOTI ( J-1 ) =0.0 
1=N 




370 


A ( 1 ) =A ( 1 ) * X1«A ( 1 «1 ) 
1 = 1-1 

1F(1-J)380.370,370 




380 


J=J«1 




15 


CONTINUE 

no fio j=j,n,2 

1TRV=0 

IF(J.LT.N)GO TO 41 
P=A (N) 

0 = A(N-n 
GO TO SO 




41 


1F(N0R0.GT.0)G0 to 43 




C 


TRY PREVIOUS ROOT FOR STARTING 
P=PJ»R00TR ( J) 


POINT 




0=0J* (ROOTR( J) *ROOTR( J) .ROOTI ( J) *ROOTI ( J) ) 
GO TO 141 


43 


P = 0.0 
0=1TRY 

IFIITRY.GT.DGO to 47 
ITRY=ITRY, 1 




1^1 


DO 451 ITER=1 , ITMAX 

80=0.0 

B1=0.0 

C0=0.0 

C1=0.0 

C2=0.0 




c 


— COMPUTE COLUMN OF 8*S AND C»S 
I=N»1 




20 


00=A ( 1 ) -PeBl-Q*B0 
B0=B1 
81=00 
1 = 1-1 

IF(I-J»2)40,40,30 




30 


00=B1-P«C2-0*C1 
C0 = C1 
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Cl*C2 
C3=00 
GO TO 20 



C SPECIAL roRM FOR LAST COEFFICIENT 

40 C2=-P*C1-0»CO 



0D=C1*C1-C2*C0 

OP=PO*C1"81*CO 

00=R1*C1-P0*C2 

IF ( a3S (DO) ,LE. ABS(T0L2«DP> .OR. ABS(OD) .LE. ABS (T0L2«00) )GO TO 43 
0P=0P/00 
00 = 00/00 
P=P.DP 
0 = 0.00 

ERR=DP*OP*DO»DO 

CR1T=P«P»0«0 

ir(CRIT.LE.TOL2*ERR)GO TO 43 
CRI T=ERR/CRI T 



ie(Iter.gt.t.and.crit.lt.tol3)go to 50 

C 00 AT least ITMIN iterations 

IF(ITEP-ITmiN)45.45.44 
C IF ERROR IS decreasing. KEEP GOING 

44 IF(CRIT.GE.CRIT0)G0 to 46 
IF (crit.lt.tolt)Go to so 

45 eRRO=ERR 
CRITO=CRIT 

451 CONTINUE 

C DROP OUT rfHEN ERR STARTS TO INCREASE OR ITER IS GREATER THAN ITMAX 

46 IF (CRITO-TOL) 50.50.43 

C IF ACCURACY IS insufficient, GIVE ERROR OUTPUT AND PROCE'^O 

47 IR=3 

WRITE (5.5000) EPRO.ERR. CRI TO. CRI T 
SOOO FORMAT ( IP. 10G12.4) 

SO PP=-P/2,0 

OISC=PP®PP-Q 

IF(DISC)52,54.S4 
52 DISC=SORT (-DISC) 



ROOTR( J-I ) rPRonOM 
ROOTR( J) =pooDOM 
ROOTI ( J-1 ) =-OISC»OOM 
ROOTI ( J) =DISC«DOM 
GO TO 56 

54 OISC=SOPT (DISC) 

ROOTR( J-I ) = (PP-DISC) *00M 
ROOTR( J) - (PP.DISC) «DOM 
ROOTI (J-I ) =0.0 
ROOTI (J) =0.0 

C synthetic DIVISION WITH QUADRATIC FACTORS 

56 I=N 
JL = J 

A (N) =A (N) -P 
60 1=1-1 

IF(I-J)80.70.70 

70 A(I)=A(I)-R*A(I.1)-0«A(I*2) 

GO TO 60 
80 CONTINUE 

IF(JL-N) 90.100. loo 
90 ROOTR (N) =-A (N) 

ROOTI (N) =0.0 
100 CONTINUE 
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200 IER=IR 
RETURN 
END 
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APPENDIX C - PROGRAM USER'S MANUAL 



This appendix will discuss the calculations performed in 
the various subroutines, the output format options, and the 
input data deck. 

DESCRIPTION OF PROGRAM 

main - The Main program does some simple calculations, 
contains formats for I/O, determines from the I/O options 
specified what calculations are to be performed, and calls 
the appropriate subroutines to perform the calculations. 

SUBROUTINE UGLYDK - This subroutine takes the five 
inputed offsets and uses a spline-cubic curve fitting routine 
to interpolate ten offsets in equal increments of T'/IO from 
the base line to T' . The subroutine is called once for each 
set of inputed offsets. The subroutine requires solution to 
a set of simultaneous linear equations which is done by the 
IBM SSP/SIMQ subroutine which should be available to most 
users . 

SUBROUTINE SPOO - This subroutine performs the calculations 
necessary to find the form coefficient k. For values above 
k = 0.19 equation (6a) is used to determine k. This subroutine 
performs the mapping necessary to determine the regression 
analysis form coefficients. This subroutine is called once 
for each set of inputed offsets. The subroutine requires the 
roots of a fifth degree polynomial which is accomplished in 
the BAIRS subroutine. 
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SUBROUTINE SLAMPR - This subroutine takes the relative 



motions, relative velocities, form coefficients, significant 
wave height, and drafts and calculates the following: 

1. number of slams at each half-station, 

2. longitudinal extent of slam, 

3. probability of slam at each half-station, 

4. most probable extreme pressure at each half-station, 

5. extreme pressure for design at each half-station, 

6. and the most probable and extreme velocities of 
impact at each half-station. 

This subroutine uses the various equations developed in Section 
IV to determine the above information. 

SUBROUTINE RELMV - This subroutine takes the non-dimen- 
sional relative motions and relative velocities inputed at 
station 1 and 4 from the tables and uses linear interpolation 
to find the values for the 17 half-stations. Seventeen half- 
stations were considered since slam damage has been reported 
as far as 40% aft of the forward perpendicular. Linear inter- 
polation is valid within the linear theory used in seakeeping 
programs . 

SUBROUTINE FORCES - This subroutine performs the temporal 
integration of pressure determined by SLAMPR and PRESS sub- 
routines. The output is the pressure time history for each 
half-station in the slam region. 
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SUBROUTINE PRESS - This subroutine takes the pressure at 
the keel from SLAMPR and calculates the girth distribution of 
vertical pressure assuming linear variation from the flat to 
T' and resolves the vertical component at each elevation using 
the hull offsets to calculate the angle 3 that the normal 
makes with the vertical (Figures 8a & 8b) . 

SUBROUTINE KEFFEC - This subroutine is not used in any 
of the calculations but determines an average hull form coeff- 
icient. This value is a measure of a ship's overall suscep- 
tibility to slam forces. For two similar ships the one with 
the larger average k will have larger slam forces and impulses. 
The average is done over the forward eleven half-stations 
since slamming rarely occurs aft of this point. 

SUBROUTINE APPROX - This subroutine uses Simpson's 
integration to find the slam impulse by integrating the 
pressure time history from FORCES at each half -station . This 
subroutine calls the IBM SSP/QSF subroutine to perform the 
integrations. Once the impulse is known, the amplitude of 
an equivalent half-sine pulse with equal duration can be 
found. This sine approximation is determined to input into 
the S.H.V.R.S. dynamic analysis program. 

SUBROUTINE MAXFOR - This subroutine fits a parabola 
through the three highest force points at each half-station and 
the first derivative is set equal to zero to find the maximum 
force at each half-station. 
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SUBROUTINE PL OT - This routine uses the M.I.T. Joint 



Computer Facility PICTR plotting routine to plot the force 
time history for each half-station on one graph. Figures 9 
to 12 are computer plots using this routine. This is the only 
subroutine that would require changing if running the program 
on another computer. 

OUTPUT OPTIONS 

There are several output options available from the 
program depending on the information desired. Options are 
specified in the data deck by assigning integer values to 
integer variables. The following options are available: 

1. 101=1 This option prints information from the SPOO 
subroutine to aid in determining problems in this 
routine. See Figure 22 for the information printed. 
101 = 0 Delete the above output. This is the 
normal option. 

2. 103=1 This option allows a plot of the pressure 
time history. See Figures 13 to 16 for typical plots. 

103 = 0 Delete the plots. 

3. 104 = 1 This outputs a table of inputed offsets and 
the ten calculated offsets for each half-station. It 
is desirable to specify this option the first time a 
set of lines are read so as to check the data. See 
Figure 23 for an example output. 

104 = 0 Delete table of offsets. 
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- 1 This option specifies that the calculations 
and output values will use most probable extreme 
pressure values. 

106=0 This option specifies that the program uses 
the extreme pressure for design. 

T 07 = 1 The force time history is outputed in 
tabular format. See Figure 24 for an example of 
the output information. 

107=0 Delete table of forces. 

6. 1010 = 0 Write table of information for each half- 
station in the slam region. See Figure 25 for the 
information available in this option. 

1010 = 1 Write a single table of information for 
all the stations in the slam region. Tables 7 to 10 
were prepared using this option. 

1010 = 2 Delete both of the above tables. 

Any combination of the above options may be specified to get 
the information desired. Since the main cost of running the 
program is the output cost, it was felt that being able to 
eliminate unnecessary output information was necessary. 
Several input options are available: 

1. 105 = 1 This allows the inputing of all the relative 
motions and velocities from seakeeping results if 
available . 

105 = 0 Input relative motions and velocities for 
stations 1 & 4 only. (Non-dimensional values from 
Tables [15] ) . 
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2. 



Time of storm may be specified if desired. If 
TSTORM is specified, it will be used in the calcul- 
ations. If set equal to zero, the TSTORM values will 
be obtained from Figure 4. 

3. The traveling velocity may be specified if desired. 
The value specified should be for a 520 foot ship 
and the program will Froude scale this figure. If 
the traveling velocity (VAFT) is set equal to zero, 
the traveling velocity will be calculated by the pro- 
gram as discussed in Section IV. 

4. The Froude number or the speed in knots can be 
specified. Both need not be inputed. 

DATA DECK 

In this section the input format for the data deck will 
be discussed. 

CARD SET #1 

CARD 1 - Format (3 I 5) - 101, 102, 104 

The values for lOl and 104 are discussed above. 

102 is always set equal to 1 and indicates that 
ships lines follow. 

CARD 2 — Format (20 A 4) — Title and information on ship. 
This card contains any information desired by the 
user to designate the ship being run. This data is 
printed in the output as a heading only and is not 
used by the program. 
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CARD 3 - Format (5 F 10.5) - BB, HH, ALBP, ALWL, STASP 
BB = Half Beam 
HH = Design Draft 

ALBP = Length Between Perpendiculars 

ALWL = Length of Load Waterline 

STASP = Station Spacing (not half-stations) 

CARD SET #2 

This card set contains the ships lines necessary for 
slam calculations. 

CARD 1 - Format (5 F 10.5) 

This card contains the waterlines at which the 
offsets were taken. The first value must be 
the 0.0 foot waterline. 

CARD 2 - Format (5 F 10.5) 

The offsets for the five waterlines in the same 
order as CARD 1 for the first half-station. 

CARD 2 is repeated for as many half-stations as 
desired up to a total of 17 half-stations. 

CARD 3 - BLANK CARD 

This card indicates that the last set of offsets 
has been inputed. 

CARD SET #3 

This set contains only one card. If no plots are desired 
(103 = 0 ) , this card is left blank. 
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Format (I 5, 4 F 10.5) - IPLT, XSCL 

IPLT = +1 - This specifies that the plot routine 
scales the plots. (AUTOSCALE) 

IPLT = -2 - This specifies to the plot routine that 
the values in XSCL be used to scale the plots. 
XSCL(l) = 0.0 = Time slam starts 
XSCL (2) = End of time axis 
XSCL (3) = 0.0 

XSCL (4) = Maximum force on vertical axis. 

If IPLT = +1, XSCL need not be specified. This same 
data will be used to scale all plots until a new set of lines 
are read. 



CARD SET #4 

Card set #4 contains the data necessary for a particular 
set of operating conditions. This set is repeated for as many 
operating conditions as desired. 

CARD 1 - Format (6 I 5) - 109, 103, 105, 106, 107, 1010 

109 = 1 always = 1 - indicates a new set of operating 
conditions follow. The other integer variables 
were discussed above. 

CARD 2 - Format (20 A 4) 

This card contains information which is printed on 
the output to distinguish this operating condition 
from the others. 
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CARD 3 Format (8 F 10.5) - FROUDE, SPEED, SEA, HF, HA, 

CB, VAFT, VTHR 

FROUDE = Froude number 

SPEED = Ship speed in knots 

SEA = Significant wave height (feet) 

HF = Draft Forward (feet) 

HA = Draft Aft (feet) 

CB = Block coefficient of underwater portion of 
hull for this draft. 

VAFT = Traveling velocity for 520 foot ship (ft/sec) 
VTHR = Threshold velocity for 520 foot ship (ft/sec) 
CARD 4 - (F 10.5) - TSTORM 

TSTORM = Duration of storm in hours 
CARD 5 - IF 105 = 0 , Format (4 F 10.5)- RM(2), RV(2), 

RM ( 8 ) , RV ( 8 ) 

RM(2) = relative motion from tables at station 1 
RV(2) = relative velocity from tables at station 1 
RM(8) = relative motion from tables at station 4 
RV(8) = relative velocity from tables at station 4 

IF 105 = 1 Format (8 F 10.5) for each card - use 
as many cards as necessary to specify the non-dimen- 
sional relative motions and velocities for as many 
half-stations as offsets read in Card Set #1. The 
values are specified in pairs beginning at half- 
station 1 (i.e., RM(1), RV(1), RM(2), RV(2), ) 
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CARD 6 - IF 103 = 0 — Delete Card 

IF 103 = 1 -- FIRST 40 spaces contain X-axis 
label 

NEXT 40 spaces contain Y-axis 
label 

This card is used by plotting subroutine. 

At the end of the operating conditions for this ship a 
blank card is inserted to indicate the end of the calculations 
for this set of lines. The above card sets 1 to 4 can be 
repeated for as many ships as desired. After the last ship 
has been read, another blank card is inserted to indicate the 
end of the data deck. Therefore the last two cards in the deck 
will be blank cards indicating the end of the operating condi- 
tions and that the last ship has been read. 
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FIOURE 22- INTERMEDIATE OUTPUT FROM SPOO SUBROUTINE ( I01 r 1) 
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FIGURE 23 - LIST OF IMPUTED AND CALCULATED OFFSETS (104=1) 
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FIGURE 24 - FORCE TIME HISTORY (107=1) 
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********************************»t******* 



* HAI F STAT ion <t * 

*■ * 
^,**********■km*^^****^^**M***m**Am*■k*»*^imm^i 

rXTPEMF PPESSiiPE. = ?AA. PSI 

•-.OST PPORARLF DPrRSLiPF = 16^-. PSI 

rOPX C"IEF. (K) = pi.(A743A 

vOM-ni-^-EMSlONAL RELATIVE MOTION = Pi. 71238 
vCN-OI^'PnjRIONAL relative VELOCITY * M.08FI7I 
^.nvurp OP RLA^S O'l'^ING EXTPE'-''F STrF-''-^ = A788. 

oPO-^AFMLlTY OF SLAM IMPACT = Pl.?7F^A? 

cxT^EMP RELATIVE VELOCITY = 57*3 FI/SFC 

V.OCT PW’OPARI.F RELATIVE VELECITY = A7.2 FT/StC 

mAVIM'J'-' value 0^ Si. AM FORCE e 3543. PI WIPS 

^Mpi|(_cr VALUE = 183. KTF-SFr 

aMPLIT'ipF op simp APPPOXIMATTOM = 2372.7: KIPS 

r,!.'»ATI''N 0*^ Siam ]mPACT FoPCE = F. 12/92 SEC 
tJmt SIAM STARTS = <5'. 14886 5^0 



FIGURE 25 -INFORMATION PRINTED AT EACH HALF-STATION 

(I01<t = <i» 
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APPENDIX D 



EXAMPLE PROGRAM 



An example program data deck and sample output are 
included to aid in using this program. MARINER case 2 is 
the example run which is the fully loaded case at 15.5 knots. 
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// XEO SLAMCW 
0 1 0 







MARINFR 








AS BUILT 




3r.O 




29.75 


528.0 


537.13 


26. 






0,00 




1.00 


2.00 


3.00 


4.00 






0.5 




1.5 


1.75 


2.05 


2.3 






1.00 




3.00 


3.5 


4.1 


4,6 






1.0 




3.45 


4.3 


5.15 


5.6 






l.o 




3.9 


5.1 


6.2 


6.6 






1.5 




4.95 


6.3 


7.5 


8.05 






?.6 




6.0 


7.5 


B.6 


9.5 






4.0 




7.5 


9.2 


10.6 


11.6 






5.4 




9.0 


10.9 


12.4 


13.5 






6.95 




11.0 


13.1 


14.8 


14.0 






8.00 




12.00 


15.00 


17. 


18.2 






8.58 




15. 


17.6 


19.7 


21.0 






9.17 




18.00 


20.2 


22.3 


24.08 






13.6 




21. 


23.10 


24.9 


26.5 




blank 




0,0 


0.5 


0. 


0 5000 


. 






1 


1 


0 0 


0 


1 








2 


MARINEP - F=0. 


20 S=0. 


047 C3=0.613 L/B= 


6.84 B/T=2.8i Fully 


LOADED 






15.5 


0.047 


27.0 


27.0 


0.613 0.0 


12.0 


0.0 
















0.0280 


0.095 


0.0191 


0.069 








2 


mariner 


TIME 


(9FC) 




FORCE (KIP5) 





blank 

blank 



// END 



FIGURE 26a > SAMPLE PROGRAM ■ DATA DECK 
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this 

ON 



PROTiRAM CALCULAirc THE 

THf LOWER 1/10 or The e 



SLAM roRCE distribution vs time 
HIP from station 8 Forward. 



uAL^ e^AM 


• 38 


.00r00 


DESIGN draft ■ ?9. 75000 




• 

Q 

CO 


S88 .00 




LFN3TH LOAD WATER LINE ■ 


537.13 


5I«!TANCE 


between 


stations - 


?/>.00 





(ALL distances in FEET) 



average hull shape factor « 0.11909 



FIGURE 26b 



SAMPLE PROGRAM 



OUTPUT 
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FIGURE 26b - SAMPLE PROGRAM - OUTPUT (cont.) 
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FIGURE 26b- SAMPLE PROGRAM - OUTPUT (cont.) 




10 

6 



m 

6 



ni 

<5 



d 



o 

d 



i 

I 



m 



127 



FIGURE 26b - SAMPLE PROGRAM - OUTPUT (cont.) 





# 




I 



4 ., 



1 



1 



[ 



i 



•1 






I 

I 



f 



Thesis 16 d003 

J22223 Walker 

Slam loading for 
use in hul 1 gi rder: 
response analysis. 
98JUM*rfi •"^DISPLAY/ • 




Thes i s 165003 
V/22223 Walker 

Slam loading for 
use I n hul 1 gi rder: 
response analysis. 



thesW22223 

hull girder 




3 2768 001 92880 7 

DUDLEY KNOX LIBRARY 



